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Introduction 


Background 


Unexploded  Ordnance  (UXO)  cleanup  is  the  number  one  priority  Army 
Cleanup  Problem  requirement  and  is  identified  as  a  major  problem  throughout 
the  Department  of  Defense  (DoD).  The  problem  is  enormous  in  scope  and 
complexity.  In  excess  of  1 1  million  acres  (approximately  44,500  km^)  of  land 
are  contaminated  with  UXO.  The  UXO  ranges  in  size  from  20-mm  rounds  to 
2,000-lb  bombs.  Approximately  50  percent  of  the  UXO  is  estimated  to  lie  on  or 
very  near  the  surface,  while  the  remaining  50  percent  is  buried  at  depths  to  6  m 
or  more.  Before  the  UXO-contaminated  DoD  lands  can  be  cleaned  up 
(remediated),  the  extent  of  UXO  contamination  must  be  determined;  i.e.,  the 
surface  and  buried  UXO  must  be  located.  To  achieve  the  UXO  location 
requirement,  Johnson  et  al.  (1996)  define  a  hierarchical  four-stage  systems 
process  of  UXO  sensing,  consisting  of  (1)  prescreening,  (2)  cuing,  (3)  detection, 
and  (4)  classification.  Prescreening  is  the  stage  of  prioritizing  potentially  UXO- 
contaminated  lands  for  follow-on  survey  based  on  historic  records  and 
subsequent  intended  land  use.  Of  the  lands  prioritized  for  follow-on  survey,  the 
cuing  stage  consists  of  location  of  regions  of  UXO  fields,  within  the  larger 
areas,  by  an  expedient  large-area  coverage  approach  (likely  an  airborne  survey). 
It  seems  logical  to  combine  the  prescreening  and  cuing  stages  of  Johnson  et  al. 
(1996)  into  a  single  screening  stage.  In  this  report,  an  additional  stage  of 
discrimination  is  defined.  The  term  identification  is  also  defined  and  is 
considered  the  same  as  the  classification  stage  of  Johnson  et  al.  Thus,  following 
the  example  of  Johnson  et  al,  this  report  considers  a  four  stage  UXO  location 
process:  screening,  detection,  discrimination,  identification/classification. 

While  surface  UXO  can  be  located  visually  and  with  remote  (airborne)  imag¬ 
ing  techniques  in  many  cases,  the  surface  area  to  be  surveyed  is  very  large  and 
cluttered  with  metallic  and  oAer  cultural  debris.  In  addition,  the  vegetative 
covers  of  the  UXO  contaminated  lands  vary  greatly  as  does  the  surface  soil  and 
rock  type  and  texture.  Although  the  surface  UXO  problem  is  large  and 
complex,  the  fact  that  it  can  be  addressed  with  remote  imaging  systems  is  a 
major  advantage,  since  such  systems  can  survey  large  areas  rapMy  with  high 
resolution.  Demonstrations  of  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  Remote  Minefields  Detection  System  (REMIDS)  airborne  system 
have  successfully  located  and  mapped  surface  UXO  (&nnett  1995).  Since 
surface  UXO  and  large  concentrations  of  surface  ordnance  debris  are  indicators 
of  the  possible  presence  of  buried  UXO,  airborne  imaging  survey  systems  such 
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as  REMIDS,  possibly  augmented  with  additional  sensors,  can  make  a  major 
contribution  to  the  screening  stage  for  location  of  buried  UXO  fields. 

The  second  stage  in  the  sensing  hierarchy  involves  the  detection  of  buried 
UXO.  General  purpose  detection  of  buried  UXO  requires  the  application  of 
ground  surface  geophysical  surveys  to  detect  geophysical  anomalies  caused  by 
the  buried  UXO.  Detection  of  buried  UXO  is  identified  as  a  major  technology 
shortfall  in  numerous  recent  studies  and  field  demonstrations  (e.g.,  Johnson 
et  al.  1996,  U.S.  Army  Environmental  Center  (USAEC)  1994,  Office  of  the 
Secretary  of  Defense  (OSD)  1997),  and  is  a  considerably  more  difficult  problem 
than  surface  UXO  mapping.  Discrimination  of  anomalies  likely  caused  by 
ordnance-like  targets  from  “false  alarm”  anomalies  caused  by  other  buried 
objects  (particularly  metallic  objects)  and  geologic  features  is  even  more 
difficult. 

Past  research,  development  and  demonstration  efforts  to  address  the  buried 
UXO  problem  have  been  disjointed  and  have  concentrated  on  adaptations  of  off- 
the-shelf  technology  and  demonstrations  at  test  sites.  The  present  work  is  the 
initial  investment  by  the  Strategic  Environment  Research  and  Development 
Program  (SERDP)  to  address  technology  shortfalls  in  buried  UXO  detection  and 
discrimination.  On-going  and  new  SERDP  efforts  will  continue  the  investment 
in  detection  and  discrimination  and  initiate  study  of  UXO  identification.  UXO 
identification  (classification)  is  the  final  step  in  UXO  sensing,  and  consists  (a)  of 
determining  the  specific  type  of  ordnance  that  most  likely  produces  a  given 
ordnance-like  anomaly  and  ultimately  (b)  of  confirming  the  presence  of 
explosives  associated  with  buried  ordnance-like  objects.  This  report  documents 
initial  efforts  to  develop  procedures  for  data  integration  and  interpretation  of 
multisensor  (multimethod)  datasets  for  UXO  detection  and  discrimination. 


UXO  Detection,  Discrimination,  and 
Identification:  Status 


The  status  of  capability  for  buried  UXO  detection,  discrimination,  and  iden¬ 
tification  is  summarized  as  follows:  (a)  can  detea  UXO,  within  definable  limits; 
0)  cannot  effeaively  discriminate  UXO  anomalies fi-om  “false  alarm” 
anomalies;  (c)  cannot  identify  UXO.  The  definable  limits  for  item  a  refer  to 
combinations  of  ordnance  size  and  burial  depth  that  result  in  geophysical 
anomalies  at  the  surface  which  can  likely  be  detected  relative  to  site-specific 
background  noise  (geologic  background  and  cultural  clutter).  “False  alarm” 
anomalies  are  caused  by  buried  ordnance  debris,  other  metallic  objects,  gravel 
and  cobbles,  soil  heterogeneities,  tree  roots,  and  other  natural  and  cultural  fea¬ 
tures.  Without  significant  discrimination  capability,  large  numbers  of  fdse 
alarms  that  must  be  verified  (dug  up)  are  the  dominant  cost  and  time  drivers  for 
UXO  site  cleanup  (remediation). 

The  results  of  recent  field  demonstrations,  such  as  the  Jefferson  Proving 
Ground  (JPG)  Technology  Demonstrations  (USAEC  1994,  1995,  1996,  1997; 
Altshuler  et  al.  1995),  exhibit  buried  ordnance  detection  probabilities  exceeding 
90  percent  by  Phase  III  (as  shown  in  Tables  1  and  2).  However,  even  with 
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Table  1 

A  Synopsis  of  JPG  UXO  TD  Phase  II  Results 

(160  UXO  Targets  Buried  for  Phase  Hi 

Demonstrator 

Targets 

Reported 

Ordnance 

Detection 

Rate  (%) 

*False  Alarms*' 

Per  Hectare 

No.  False  Alarms 

Per  Ordnance 

Item  Detection 

Geophex 

398 

71 

19.7 

3.41 

Geometries 

521 

83 

26.9 

3.96 

Parsons 

602 

85 

32.5 

4.68 

Bristol 

566 

62 

38.3 

6.97 

ADI  (Combined) 

598 

65 

34.5 

9.35 

Coleman 

280 

29 

15.9 

9.56 

Scintrex 

255 

50 

45.3 

10.10 

GeoPotentia 

168 

11 

12.0 

13.00 

Geo-Centers 

1,409 

72 

84.0 

20.70 

Vallon 

1,903 

57 

225.9 

68.00 

Table  2 

A  Synopsis  of  JPG  UXO  TD  Phase  III  Results 

Artillery  and  Mortar  Range  (Scenario  2) 

(117  Targets  -  67  Ordnance;  50  Nonordnance) 

Demonstrator 

Targets 

Reported 

Ordnance 

Detection 

Rate  (%) 

"False  Alarms" 

Per  Hectare 

No.  False  Alarms 

Per  Ordnance 

Item  Detection 

NAEVA 

202 

97 

19.0 

1.37 

174 

67 

21.1 

2.20 

Geometries 

282 

90 

38.4 

3.00 

Ensco 

279 

70 

43.6 

4.34 

Geo-Centers 

486 

93 

80.7 

6.10 

ADI 

456 

85 

76.8 

6.32 

Rockwell 

151 

21 

27.1 

9.07 

Geo  Potential 

23 

3 

4.3 

10.00 

GRI 

1,319 

90 

258.2 

20.15 

ordnance  detection  improving  to  acceptable  rates,  the  number  of  false  alarms  is 
unacceptably  high,  i.e.,  poor  discrimination  capability.  For  JPG  Phase  H,  four 
demonstrators  had  ordnance  detection  rates  >  70  percent;  the  number  of  felse 
alarms  for  each  ordnance  item  detected,  however,  ranged  from  3.4  to  20.7  for 
these  demonstrators  (Table  1).  Although  JPG  Phase  III  was  considerably  easier 
for  ordnance  detection  than  Phase  n  in  that  the  ordnance  items  were  consistently 
shallower  (Figure  1),  it  is  notable  that  ordnance  detection  rates  improved 
considerably.  Four  demonstrators  for  JPG  Phase  HI  Scenario  2 
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Figure  1 .  Comparison  of  ordnance  target  burial  depths  for  PG  Phases  II  and  III 


(Table  2)  had  ordnance  deteaion  rates  i  90  percent;  but  the  numbers  of  false 
alarms  per  ordnance  item  detected  ranged  from  1.4  to  20.2,  still  unacceptably 
high  although  showing  some  improvement.  The  JPG  and  other  field  demonstra¬ 
tions  exhibit  limited  capability  for  ordnance  identification  or  classification  (item 
c  above).  Even  classification  into  broad  ordnance  categories,  such  as  bombs, 
projectiles,  and  mortars,  is  not  reliably  possible  with  current  fielded  detection 
systems.  Capability  for  verification  of  explosive  content  in  buried  ordnance 
does  not  exist. 


UXO  Detection,  Discrimination,  and 
Identification:  Approaches 


Detection  of  UXO  can  be  achieved  with  single  or  multiple  geophysical 
method  (single-  or  multisensor)  approaches.  Single-sensor  (total  field  magne¬ 
tometers  or  electromagnetic  induction  systems)  approaches  have  demonstrated 
probabilities  of  ordnance  detection  >  75  percent,  however  most  demonstrators 
use  a  multisensor  approach.  Generally,  demonstrators  with  considerable  field 
experience  with  geophysical  methods  and  the  interpretation  of  geophysical  data 
exhibit  consistently  better  UXO  detection  performance  in  both  demonstrations 
and  live  site  surveys.  Geophysicists  have  long  recognized  the  value  of  inte¬ 
grated,  complementary,  multimethod  approaches  for  a  variety  of  subsurface 
detection  and  mapping  objectives.  Complementary  geophysical  methods  mea¬ 
sure  parameters  affected  by  different  physical  properties  of  the  subsurface. 
Detection  of  UXO  requires  high-quality  data  collection  efforts  with  special 
attention  to  accurate  positioning  of  measurements.  The  next  requirements  for 
UXO  detection  consist  of  a  procedure  for  selection  of  anomali^  and  a  decision 
criteria  for  target  declaration.  For  multisensor  datasets,  the  anomaly  selection 
criterion  may  be  a  requirement  for  coincident  anomalies  on  two  (or  more)  com¬ 
plementary  datasets.  The  final  requirement  is  to  interpret  the  single-  or 
multisensor  datasets  for  information  about  the  target,  e.g.,  map  position, 
estimated  depth,  and  ferrous  mass. 

There  is  a  general  consensus  that  measurements  with  presently  fielded  single¬ 
sensor  geophysical  systems  do  not  contain  sufficient  information  to  allow 
discrimination  of  UXO  anomalies  from  false-alarm  anomalies.  General-purpose 
anomaly  discrimination  requires  multisensor  data  acquisition  using  two  or  more 
complementary  geophysical  methods  or  the  use  of  emerging  geophysical 
systems  which  measure  much  more  information  about  anomalies  in  terms  of 
time  variation  effects,  frequency  variation  effects,  or  spatial  gradients  of  fields. 
The  emerging  geophysical  systems  can  be  described  as  multichannel, 
multifrequency,  and/or  multicomponent.  Indeed  a  multicomponent  measure¬ 
ment  system  can  be  considered  a  multisensor  system,  even  if  the  same  physical 
parameter  is  measured  along  different  spatial  directions.  The  most  important 
step  for  UXO  discrimination  is  multisensor  data  integration  and  interpretation 
(sometimes  called  “data  fusion”).  A  key  requirement  for  successful  multisensor 
discrimination  is  the  development  of  innovative  discrimination  algorithms.  The 
different  approaches  to  and  levels  of  sophistication  of  data  integration  and 
interpretation  are  discussed  in  Chapter  2. 

Identification  or  classification  of  UXO  anomalies  as  specific  ordnance  types 
requires  the  development  of  new  and  innovative  geophysical  technology.  A 
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possible  approach  may  be  the  detection  of  characteristic  acoustic  or  electromag¬ 
netic  resonance  responses  that  allow  identification;  this  approach  will  most 
likely  be  applicable  to  localized  interrogation  of  the  subsurface  beneath  previ¬ 
ously  detected  anomalies.  More  challenging  will  be  the  definitive  verification 
of  explosive  content  in  the  target.  Explosive  detection  will  require  that  the 
ordnance  be  either  extremely  shallow  or  that  a  sensor  probe  be  placed  in  close 
proximity  to  a  previously  detected  target  identified  as  likely  UXO.  Possible 
sensor  approaches  for  explosives  detection  include  neutron  activation,  neutron 
backscatter,  or  soil  gas  analyses.  For  convenience,  the  acronym  UXO  is  often 
used  rather  loosely  in  this  report  to  refer  to  an  ordnance  item,  whether  inert  or 
truly  an  unexploded  ordnance  (UXO).  All  of  the  ordnance  items  at  test  sites 
referenced  in  Chapter  3  (“Data  Sources  for  the  Investigation  ”)  are  inert. 


Scope  of  Report 


Chapter  2  briefly  surveys  the  geophysical  methods  commonly  utilized  for 
UXO  detection  and  discrimination,  discusses  the  fundamental  physics  and 
phenomenology  concepts,  and  surveys  the  approaches  to  multisensor  integration 
and  interpretation.  Sources  of  data  for  this  investigation  are  described  in 
Chapter  3.  Analysis  methods  for  data  management  and  detection  enhancement 
are  discussed  in  Chapter  4,  and  a  synopsis  of  current  UXO  detection  capability 
is  presented  in  Chapter  5.  Chapter  6  describes  efforts  to  develop  fundamental 
phenomenological  modeling  capability  for  UXO.  Chapter  7  presents  the  results 
of  investigations  of  UXO  detection  and  discrimination  utilizing  emerging 
geophysical  technology,  and  conclusions  and  recommendations  are  outlined  in 
Chapter  8. 
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Geophysical  Methods  and 
Multisensor  Interpretation 
Approaches 


Geophysical  Methods  for  UXO 
Detection  and  Discrimination 


For  completeness,  a  brief  overview  follows  of  the  geophysical  methods,  con¬ 
cepts,  and  applications  rationale  for  buried  UXO  detection  and  discrimination. 
Currendy,  most  live  site  UXO  detection  surveys  utilize  total  field  magnetic 
and/or  electromagnetic  (EM)  induction  methods.  Often  the  magnetic  surveys 
include  measurement  of  the  total  field  vertical  gradient.  Demonstrations  at  test 
sites  have  included  ground  penetrating  radar  surveys,  various  emerging  electro¬ 
magnetic  induction  systems,  various  emerging  magnetic  systems,  and  airborne 
surveys.  Generally,  ground  penetrating  radar  (GPR)  surveys  for  UXO  detection 
have  been  unsuccessful  due  to  inappropriate  site  conditions  (high  EM 
attenuation  and/or  high  background  clutter)  for  application  of  radar, 
inappropriate  selection  of  radar  system  parameters,  inexperienced  personnel,  or 
a  combination  of  these  factors.  Airborne  surveys  for  detection  of  buried  UXO 
(say  >  10  cm  deep)  have  been  spectacularly  unsuccessful;  evidence  from  the 
JPG  demonstrations  failed  to  indicate  any  detection  capability  for  airborne 
systems  (Altshuler  et  al.  1995,  USAEC  1996).  Additional  geophysical  methods 
proposed  for  application  to  UXO  detection  and  discrimination  include 
microgravimetry  and  acoustic/seismic  systems. 

Geophysical  methods  interrogate  or  probe  the  physical  properties  of  subsur- 
&ce  materials  and  specifically  detect  the  presence  of  contrasts  in  physical  prop¬ 
erties  of  subsurface  materials.  The  geophysical  techniques  are  classified  in 
Figure  2  based  on  the  nature  of  the  “energy  source*  for  the  method  and  on  the 
nature  of  the  physical  phenomenon  that  is  exploited.  The  potential  fields  meth¬ 
ods,  for  example,  are  generally  passive  and  involve  measurement  of  natural  or 
preexisting  potential  fields,  e.g.,  the  earth’s  gravitational  and  magnetic  fields. 
Practical  application  of  the  potential  fields  methods  detect  anomalies  in  the 
fields  caused  by  “localized”  contrasts  in  physical  properties  of  subsurface  mate¬ 
rials  (e.g.,  density  contrast  for  the  gravity  method).  EM  methods  are  both 
active  and  passive;  although  methods  applicable  to  UXO  deteaion  are  active,  in 
that  an  EM  transmitter  is  part  of  the  system. 
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Nature  of  the  Source 


Active  Methods  ■■■ 


SEISMIC  METHODS 
ELECTRICAL  RESISTIVITY 
EM  METHODS 
RADIOACTIVITY  METHODS 


Passive  Methods 


GRAVITY  METHODS 
MAGNEVC  METHODS 
SELF  POTENTIAL  METHOD 
MAGNETOTELLURIC  METHOD 
RADIOACTIVITY  METHODS 


Nature  of  the  Physical  Phenomenon 

POTENTIAL  FIELDS  METHODS 
EM  METHODS 
ELECTRICAL  METHODS 
SEISMIC  METHCJOS 
RADIOACnVITY  METHODS 


Figure  2.  Classificadons  of  geophysical  methods 


Measurements  during  geophysical  surveys  are  generally  in  basic  units,  e.g., 
distance,  time,  voltage,  current,  magnetic  field  strength,  or  meter  units.  The 
process  of  converting  measured  quantities  to  interpreted  parameters,  in  terms  of 
a  subsurface  model,  an  image  of  the  subsurface,  or  surface  maps  of  individual 
parameters,  is  illustrated  in  Figure  3.  Proceeding  from  the  calculated  param¬ 
eters  to  the  interpretation  stage  requires  some  form  of  modeling  and  an 
understanding  of  fundamental  concepts  and  phenomenology.  While  it  may  be 
possible  to  produce  surface  anomaly  maps  direcdy  after  the  measurement  or 
calculated  parameters  stages,  it  is  not  possible  to  achieve  a  rigorous  multi¬ 
method  interpretation  without  fundamental  phenomenological  understanding  of 
the  methods.  A  synopsis  of  the  major  geophysical  methods,  concepts  (key  prin¬ 
ciple),  measured  quantities,  the  key  physical  parameters,  and  the  interpreted 
parameters  is  shown  in  Figure  4.  Finally,  the  tabulation  below  lists  the 
geophysical  methods  likely  to  be  encountered  in  UXO  detection  and  discrimi¬ 
nation  efforts  and  gives  abbreviations  which  will  be  used  subsequently  in  this 
report: 


Geophysical  Methods  and  Abbreviations 


I.  -  Electromagnetic  Induction  Methods 

FDEM  -  Frequency  domain  EMI;  quadrature  and  in-phase  component  mea¬ 
surement;  single-frequency;  multifrequency;  multicomponent 
TDEM  -  Time  domain  EMI;  transient  decay  signatures;  single-channel 
(gate);  multichannel;  multicomponent 

II.  Magnetic  Methods 

TFM  -  Total  field  magnetic 

TFG  -  Total  field  vertical  magnetic  gradient 

III.  GPR  -  Ground  Penetrating  Radar 

Pulse  (time  domain),  CW  (frequency  domain),  chirp,  etc.,  systems 

IV.  Microgravity  Methods 

V.  Seismic  or  Acoustic  Methods _ 
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MEASUREMENTS 


CORRECTIONS 
SIMPLE  RELATIONS 
CORRELATIONS 
INSTRUMENT  FACTORS 


MAGNETIC  FIELD  STRENGTH 

DISTANCES 

TIMES 

RESISTANCES  (  VOLTAGES,  CURRENTS  ) 
METER  UNITS 


CALCULATED  PARAMETERS 


MODELING 


FORWARD 
INVERSE 
EMPIRICAL  CORRELATIONS 
QUALITATIVE 


DEPTHS 
GEOMETRY 
THICKNESSES 
DENSITIES 
SUSCEPTlBILmES 
RESISTIVITIES 
coNDUcnvmEs 
SEISMIC  VELOCITIES 
ELASTIC  MODULH 


SPEEDS  (  VELOCITIES  ) 
RESISTIVITIES 
ACCELERATIONS 
coNDucnvmEs 
MAGNETIC  FIELD  STRENGTH 


INTERPRETED  PARAMETERS  - 
"Model"  of  Subsurface 


Figure  3.  Concept  of  geophysical  measurements  and  dieir  transformation  to  interpreted 
parameters 


Surveys  including  two  or  more  of  the  methods  discussed  above  are  multi¬ 
method  or  multisensor  surveys  or  approaches  to  characterization  of  the  sub¬ 
surface.  Characterization  includes  detection  and  discrimination  as  defined  for 
UXO.  The  methods  from  the  tabulation  above  (main  categories)  or  from 
Figure  4  are  termed  complementary  methods,  when  the  methods  respond  to 
different  key  physical  properties  or  sets  of  properties.  In  this  sense,  a  survey 
with  magnetic  and  EMI  methods  would  be  a  multisensor,  complementary 
methods  survey.  However,  a  survey  with  FDEM  and  TDEM  methods  would  be 
multisensor  but  not  a  complementary  methods  survey  (Butler  1986,  Butler  and 
Fitterman  1986).  Also,  a  survey  with  a  multiaxis  TDEM  system  is  considered 
multisensor  but  does  not  constitute  complementary  measurements. 
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Chapter  2  Geophysical  Methods  and 


Multisensor  Interpretation  Approaches 


Figure  4.  Geophysical  methods,  measurements,  and  interpreted  parameter; 


Multisensor  Data  Interpretation  Approaches 
for  UXO  Detection  and  Discrimination 


Classically  or  traditionally,  geophysicists  have  achieved  multisensor  (multi¬ 
method)  data  interpretation  in  an  ad  hoc,  empirical  manner,  ^ich  relied  exten¬ 
sively  on  knowledge  of  phenomenologies  and  experience.  The  data  from  each 
survey  method  are  interpreted  individually  and  then  integrated  with  each  other 
and  with  any  available  geological  ground  truth  to  form  an  interpreted  model  of 
the  subsurface  (e.g.,  Butler  et  al.  1996).  Various  types  of  analysis  and  data 
management  tools  have  been  developed  to  aid  the  interpretation  process.  The 
formal  interpretation  of  a  given  dataset  is  known  as  geophysical  inversion, 
where  an  earth  model  is  deduced  directly  from  data.  In  its  most  rigorous  form, 
geophysical  inversion  produces  a  best-fit  model  to  the  data  in  a  statistical  sense. 
Frequently  geophysical  inversion  proceeds  from  the  assumption  of  an  initial 
model;  the  physical  properties  and  dimensions  of  the  model  are  adjusted  to 
achieve  a  best-fit  to  the  data  (Parker  1994,  Meju  1994,  Butler  et  al.  1982).  It  is 
after  the  point  of  achieving  an  interpretation  of  the  individual  datasets  that 
integrated  interpretation  has  traditionally  started.  Procedures  and  capabilities 
for  achieving  true,  formal  joint  inversion  of  multisensor  datasets  are  rapidly 
emerging  (e.g.,  Sandberg  1990,  Dobroka  et  al.  1991,  Lavely  and  Grimm  1997). 
The  following  sections  survey  the  approaches  to  multisensor  data  interpretation, 
proceeding  from  the  empirical  approach  to  true  joint  inversion. 


Empirical  interpretation  procedures 
for  complementary,  multisensor  datasets 

The  ad  hoc,  empirical  approach  to  multisensor  interpretation  relies  on  con¬ 
siderable  knowledge  of  the  phenomenologies  of  the  geophysical  methods  and 
experience  of  the  interpreter.  Geophysical  survey  programs  are  conducted  with 
specific  objectives,  e.g.,  groundwater  exploration,  cavity  detection,  UXO 
detection,  etc.  The  geophysical  methods  used  for  the  programs  are  selected 
with  knowledge  of  the  physical  nature  of  the  subsurface  targets  (e.g.,  saturated 
porous  media,  air-  or  water-  or  clay-filled  cavity,  localized  metallic,  conductive 
object),  the  expected  depth  or  depth  range,  the  anticipated  target  size  and  geom¬ 
etry,  and  the  nature  of  the  geologic  material  above  and/or  surrounding  the 
target.  Data  from  each  survey  method  are  interpreted  using  well-established 
approaches  (Telford,  Geldart,  and  Sheriff  1990,  Burger  1992,  Headquarters, 

U.  S.  Army  Corps  of  Engineers  (HQUSACE)  1995).  Empirical  relationships, 
“rules-of-thumb,”  and  past  experience  with  similar  problems  then  guide  the 
multimethod  interpretation  process,  such  as  Ohm’s  law,  Archie’s  relationships, 
porous  media  density  model,  Wylie’s  time  average  equation  for  seismic  velocity 
in  porous  media,  simple  dipole  magnetic  field  model,  and  tabulations  of  seismic 
velocity,  electrical  resistivity,  magnetic  susceptibUity,  and  dielectric 
permittivity  of  geologic  materials  (see  for  example.  Sheriff  1991  and 
Carmichael  1989,  for  definitions  and  discussions  of  these  concepts). 

A  common  problem  in  geophysics  is  to  interpret  soil  and/or  rock  type  and 
condition  for  a  given  layer  in  a  layered  model  of  the  subsurface.  Seismic 
refraction  and  electrical  resistivity  survey  results  are  commonly  interpreted  in 
the  form  of  a  layered  model.  Using  the  empirical  tools  discussed  previously,  a 
qualitative  interpretation  of  the  soil  and  rock  type  and  condition  can  be  inferred 
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from  joint  consideration  of  the  seismic  compression  wave  velocity  and  the  elec¬ 
trical  resistivity.  The  basic  concepts  of  qualitative  geologic  interpretation  are 
illustrated  in  Figure  5,  although  real  geologic  conditions  may  differ  from  the 
indications  in  the  figure. 

Another  illustration  of  the  empirical  approach  to  integrated  interpretation  of 
geophysical  survey  data,  and  one  more  closely  related  to  UXO  detection  (i.e., 
localized  anomalies)  and  discrimination,  is  the  problem  of  cavity  detection 
(Butler  1994).  Subsurface  cavities  can  be  manmade  (tunnels  and  mines)  and 
natural  (caves,  caverns,  and  fissures  in  karst  areas).  The  cavities  can  be  air-, 
clay-,  or  water-filled,  and  tunnels  can  be  lined  and  contain  metal.  Common 
feamres  on  the  weathered  surface  of  limestones  in  karst  areas  are  limestone 
pinnacles  and  clay-filled  grikes  (areas  between  two  pinnacles  or  dissolution  fea¬ 
tures  in  top  of  rock;  Franklin  et  al.  1981).  Figure  6  contains  examples  of  quali¬ 
tative  interpretation  of  localized  anomalies  in  a  karst  region  geophysical  site 
characterization  (Butler  1983).  An  important  concept  of  geophysical  interpreta¬ 
tions  is  that  frequently  there  are  ambiguities  or  uncertainties.  Often  the  ambi¬ 
guity  may  be  significantly  reduced  by  consideration  of  additional  sources  of 
information.  However,  there  will  always  be  some  degree  of  fundamental  or 
inherent  ambiguity  in  geophysical  interpretations,  and  only  when  direct  investi¬ 
gation  of  the  subsurface  is  conducted,  can  ambiguity  be  totally  removed  (Butler 
et  al.  1996,  Simms  and  Butler  1992,  Simms,  Butler,  and  Powers  1995). 

Finally,  Figure  7  illustrates  the  procedure  for  qualitative  interpretation  of 
UXO-like,  localized  anomalies  in  geophysical  survey  results.  As  with  the  pre¬ 
vious  example,  there  is  ambiguity  in  Ae  interpretation  of  UXO-like  anomalies. 
Spatial  properties  of  the  anomalies  (e.g.,  spatial  wavelength)  are  used  to  esti¬ 
mate  depths  and  infer  geometries  and  orientations.  With  experience,  the  geo¬ 
physicist  (or  data  analyst)  can  become  proficient  at  discriminating  UXO 
anomalies  from  false-alarm  anomalies,  but  it  is  a  slow  and  painstaking  process. 
And  even  with  an  experienced  interpreter,  ambiguity  (false  alarms)  cannot  be 
eliminated. 


Analysis  and  data  management  tools  for 
multisensor  data  integration  and  interpretation 

Analysis  and  data  management  tools  are  key  links  between  field  acquisition 
of  large  datasets  and  integrated  multisensor  interpretation.  Approaches  vary 
from  custom-designed  analysis  and  data  management  software,  developed  to 
support  specific  sensor  systems  or  platforms  (e.g.,  McDonald  and  Robertson 
1996),  to  the  use  of  commercially  available  software,  where  the  field  data  are 
reformatted  as  necessitated  by  software  input  requirements.  Modem  graphics 
software  are  invaluable  for  visualizing  geophysical  survey  data  over  areas.  Line 
contour  plots/maps,  shades  of  gray  maps,  full  color  maps,  and  shaded  relief 
maps  greatly  enhance  capability  for  detecting  anomalies  and  data  trends. 
Multisensor  datasets,  particularly  when  collected  separately  or  at  different 
times,  may  not  be  exactly  spatially  co-located.  For  example,  EMI  and  magnetic 
survey  datasets,  even  if  collected  at  the  same  survey  line  spacing,  may  have 
different  measurement  spacings  along  the  lines.  Also,  the  survey  lines 
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Example  of  Empirical  Integration 
of  Complementary  Geophysical  Datasets 
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Example  of  Empirical  Integration 
of  Complementary  Geophysical  Datasets 

Simplistic  Assessment  UXO-Like,  Localized  Anoma! 
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Figure  7.  Empirical  integration  of  complementary  geophysical  datasets  for  UXO-like,  localized  anomalies 


themselves  may  deviate  somewhat  between  sensor  types.  Many  types  of  inte¬ 
grated  multisensor  interpretation  procedures  work  best  or  most  efficiently  with 
co-located  multisensor  data.  A  common  procedure  for  achieving  spatially 
co-located  datasets  is  to  fit  the  original  datasets  with  individual  “best-fitting” 
surfaces  and  then  sample  the  surfaces  at  regularly  spaced  grid  points.  Tbe 
subject  of  gridding  is  beyond  the  scope  of  this  report;  however,  it  is  important 
to  carefully  consider  the  gridding  procedure  and  its  possible  effects  on  trends 
and  anomalies  in  the  datasets  (e.g.,  Gallagher  1989,  Scollar  et  al.  1990, 
MacLeod  and  Dobush  1990,  Bhattacharyya  1969).  Modem  graphics  software 
can  readily  process  large  data  volumes,  fit  the  data  with  user-selectable  surfaces 
(inverse  distance  weighting,  minimum  curvature,  kriging,  polynomial 
regression,  etc.),  and  produce  a  regular  grid  of  “data”  values  from  the  surface 
fits  (e.g..  Golden  Software  1995,  Geosoft  1996).  In  this  manner,  maps  of  the 
survey  results  for  each  sensor  survey  can  be  produced  from  gridded  values  that 
are  co-located.  The  co-located  multisensor  datasets  are  then  ready  for  inte¬ 
grated  processing  and  interpretation.  Often  each  sensor  dataset  may  be  individ¬ 
ually  processed  in  some  way  (e.g.,  background  subtraction,  image  processing, 
various  kinds  of  filtering  operations,  etc.)  prior  to  integrated  processing  and 
interpretation.  Some  of  the  graphics  software  packages  (e.g.,  Geosoft  1996) 
have  relational  links  between  the  different  sensor  datasets  and  allow  direct  tran¬ 
sition  to  simple  types  of  integrated  multisensor  dataset  analyses. 

Geographic  information  systems  (GISs)  are  ideal  tools  for  management  of 
multisensor  datasets  and  all  additional  site  information  (e.g.,  topography,  vege¬ 
tation,  surface  clutter,  etc.).  Each  sensor  dataset  and  each  additional  type  of  site 
information  are  entered  as  separate  map-planes  in  the  GIS.  Importantly,  all  of 
the  map-planes  are  georeferenced.  Most  GIS’s  allow  for  input  of  any  type  of 
information  with  corresponding  position  tags  (locations).  Also,  most  GIS’s 
allow  for  the  procedures  discussed  in  the  preceding  paragraph,  i.e.,  can  fit  the 
sensor  datasets  with  best-fitting  surfaces  and  sample  at  regular  grid  points,  to 
produce  co-located  multisensor  datasets.  A  GIS  includes  relational  database 
capabilities,  so  that  one  or  more  georeferenced  map-planes  can  be  queried  to 
produce  derived  map-planes.  For  example,  a  typical  query  relevant  to  UXO 
detection  and  discrimination  applications  might  include  a  derived  map  showing 
all  locations  where  a  geophysical  anomaly  (e.g.,  magnetic  or  EMI)  is  coincident 
in  location  with  a  mapped  surface  metallic  object.  Another  example  of  a  typical 
query  is  to  produce  a  derived  map  showing  locations  (or  areas)  where 
coincident  geophysical  anomalies  occur  on  two  or  more  map  planes  (e.g., 
coincident  magnetic  and  EMI  anomalies).  Often  a  derived  map,  such  as 
described  in  the  previous  example,  will  be  the  final  product  of  an  integrated 
multisensor  interpretation  effort.  Derived  maps  may  also  be  functional  combi¬ 
nations  of  two  or  more  map-planes  (sum,  difference,  or  more  complicated  rela¬ 
tionship).  A  GIS  may  also  have  linked  or  embedded  expert  systems  or  neural 
network  classifiers  (Burrough  1986,  Millhouse  et  al.  1996). 


Integrated  multisensor  interpretation 

A  common  approach  to  multisensor  interpretation  is  the  use  of  forward 
modeling  and  comparison  to  measurements.  For  each  geophysical  method, 
forward  modeling  depends  on  postulating  a  model  of  the  subsurface  and  then 
utilizing  a  physics-based  empirical,  analytical,  or  numerical  approach  for  com¬ 
puting  the  prediaed  response  of  the  postulated  model.  The  forward  modeling 
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can  be  manually  or  automatically  iterated  to  achieve  a  fit  or  match  of  model 
prediction  to  measured  data  (a  process  known  as  geophysical  inversion).  Often 
the  forward  modeling  software  will  compute  a  goodness  of  fit  value  (e.g.,  root- 
mean-square  (rms)  error)  to  the  measured  data,  and  automated  inversion  seeks 
to  minimize  a  goodness  of  fit  measure  (Parker  1994,  Meju  1994,  Butler  et  al. 
1982).  Integrated  interpretation  occurs  when  the  forward  modeling  procedure 
includes  i/irer-method  iteration  as  well  as  intra-method  iteration  or  constraints. 

For  the  most  common  suite  of  geophysical  methods  used  for  UXO  detection 
surveys,  TFM  and  TDEM,  the  following  tabulation  indicates  some  of  the 
parameters  involved  in  models  for  analytical  or  numerical  forward  modeling  or 
as  initial  models  for  inverse  modeling; 


Magnetic  Model 

TDEM  Model 

Geometry  (e.g.,  sphere  or  prolate 
spheroid) 

Geometry  (e.g.,  sphere,  cylinder  or  prolate 
spheroid 

Length,  Diameter 

Length,  Diameter 

Depth 

Depth 

Inclination,  Azimuth 

Inclination,  Azimuth 

Magnetic  Permeability  of  UXO 

Electrical  Conductivity  of  UXO 

Magnetic  Permeability  of  Sotl/Rock 

Electrical  Conductivity  of  Soil/Rock 

Earth's  Magnetic  Field  Magnitude 

EM  Transmitter  Dipole  Moment 

Earth's  Field  Inclination  and  Declination 

Receiver  Time  Gate  (or  Gates) 

Many  of  these  parameters  are  known  for  particular  cases  and  will  be  fixed  (con¬ 
stant)  for  the  forward  modeling  and  manual  iteration  cycles.  Geometrical 
considerations  for  these  two  models  are  illustrated  in  Figure  8,  where  the  geom¬ 
etry  for  the  magnetic  and  electromagnetic  sources  and  for  the  UXO  are  shown. 
As  discussed  previously,  the  magnetic  method  is  passive,  using  the  earth’s 
magnetic  field  as  the  source,  and  the  resulting,  induced  anomalous  magnetic 
field  (superimposed  on  the  inducing  earth’s  field)  is  static.  The  EM  methods 
are  active  and  the  resulting,  induced  anomalous  electromagnetic  field  is  time- 
varying,  with  frequency  content  (FDEM)  or  time  decay  characteristics  (TDEM) 
dependent  on  the  EM  source  characteristics.  Since  the  electrical  conductivity 
and  magnetic  permeability  of  UXO  targets  are  generally  much  larger  than 
surrounding  geologic  materials,  the  induced  background  responses  from 
surrounding  geologic  materials  are  much  smaller  than  the  superimposed, 
anomalous  responses  from  the  UXO. 

The  most  common  approach  to  integrated  multisensor  interpretation 
involves  feature  and  parameter  extraction.  Features  or  parameters  are  deter¬ 
mined  from  analyses  or  inversions  of  measurements  and  then  either  used  in 
parameter  cross  plots  or  for  sequential  input  to  other  analyses  or  inversions 
(Bell  and  Barrow  1997).  For  example,  if  a  localized  anomaly  is  detected  at 
approximately  the  same  apparent  location  with  both  TFM  and  TDEM  surveys, 
the  feature  causing  the  anomalies  is  both  magnetically  permeable  and 
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Magnetic  Surveying  --  Passive  Methods  Eiectromagnetic  Surveying  -  Active  Methods 
Anomaly  Induced  by  Earth's  Magnetic  Field  Anomaly  Induced  by  Vertical  Magnetic  Dipole  Source 


Figure  8.  Geometrical  considerations  for  UXO  magnetic  and  electromagnetic  anomaly  response  modeling 


electrically  conductive  relative  to  the  surrounding  geologic  materials  and  is 
potentially  UXO.  An  integrated  interpretation  could  consist  of  (a)  determining 
the  depth  to  the  conductive  objective  from  analysis  of  TDEM  measurements  and 
then  (b)  inputting  the  depth  as  a  fixed  or  constrained  model  parameter  during 
inversion  of  the  magnetic  data.  The  previous  form  of  integrated  interpretation 
can  be  considered  a  sequential  inversion  process  or  simplified  joint  inversion;  a 
proposed  algorithm  for  accomplishing  sequential  inversion  of  magnetic  and 
TDEM  data  is  presented  later  in  this  report  (including  an  example).  A  hypo¬ 
thetical  example  of  a  parameter  cross  plot  is  shown  in  Figure  9,  interpreted 
object  mass  versus  interpreted  object  volume,  from  an  analyses  of  coincident, 
localized  gravity  and  magnetic  anomalies,  respectively.  In  principle,  a  parame¬ 
ter  cross-plot  could  be  used  for  both  discrimination  and  identification.  For 
example,  two  points  in  Figure  9,  determined  from  inversion  (interpretation)  of 
gravity  and  magnetic  data  (asterisks),  lie  very  close  to  points  representing 
known  UXO  (solid  circles),  allowing  identification  of  the  buried  object  as  a 
specific  UXO.  Another  point  in  Figure  9,  determined  from  inversion  of  gravity 
and  magnetic  data,  lies  outside  the  region  of  mass-volume  space  corresponding 
to  known  UXO,  allowing  discrimination  of  the  buried  object  as  non-UXO. 
Examples  of  various  types  of  parameter  cross  plots  are  discussed  later  in  this 
report. 


Joint  inversion 

Joint  inversion  is  the  most  rigorous  and  self-consistent  procedure  for  achiev¬ 
ing  truly  integrated  interpretation  of  multimethod  geophysical  datasets  (Meju 
1994,  Laveley  and  Grim  1997).  Rigor  and  self-consistency  do  not  necessarily 
imply  that  joint  inversion  is  the  best  or  most  efficient  or  most  accurate  approach 
to  integrated  interpretation  of  multisensor  datasets  for  UXO  discrimination  and 
identification  in  all  cases.  However,  the  “best-fitting”  model  resulting  from 
joint  inversion  of  complementary  datasets,  e.g.,  TFM  and  TDEM  data,  will 
generally  have  considerably  less  ambiguity  than  single-method  inversion  or  the 
various  forms  of  integrated  interpretation  (e.g.,  sequential  inversion)  discussed 
previously.  The  complementary  methods  interrogate  or  “illuminate”  the  target 
in  distinctly  different  manners.  For  UXO  targets,  intelligent  estimates  of 
material  property  parameters  and  of  the  contrast  between  the  UXO  and 
surrounding  materials  are  possible.  Also,  the  properties  of  the  EMI  source  and 
the  earth’s  field  are  known.  Thus,  if  the  joint  inversion  utilizes  a  common 
geometric  model  for  both  datasets,  e.g.,  a  prolate  spheroid,  the  number  of 
parameters  that  must  be  determined  in  Ae  joint  inversion  reduces  to  two 
intrinsic  properties  (length  and  diameter)  and  three  extrinsic  properties  (depth, 
inclination,  and  azimuth).  Joint  inversion  considerably  improves  the  resolution 
and  error  of  the  model  solution  (Laveley  and  Grimm  1997),  compared  to  single¬ 
method  inversion  for  the  same  number  of  model  parameters.  Joint  inversion  has 
been  attempted,  with  varying  degrees  of  success,  for  gravity  and  magnetics, 
seismic  and  electrical  resistivity,  seismic  and  GPR,  and  electrical  resistivity  and 
TDEM.  Joint  inversion  of  TFM  and  TDEM  datasets  is  a  well-posed  inversion 
problem;  however,  forward  modeling  capability  for  TDEM  for  realistic  three- 
dimensional  (3-D)  UXO  geometries  (e.g.,  prolate  spheroid)  does  not  currently 
exist. 
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Figure  9.  Cross  plot  of  parameters  derivable  from  gravity  and  magnetic  datasets  for  UXO  discrimination  and 
identification 


3 


Data  Sources  for  the 
Investigation 


Key  to  the  success  of  efforts  to  develop  multisensor/multimethod  interpre¬ 
tation  are  data  signatures  from  known  ordnance  targets.  There  is  a  sparsity  of 
published  multisensor  data,  acquired  over  known  ordnance  targets.  Initially  this 
investigation  planned  to  make  extensive  use  of  data  from  the  JPG  UXO 
Technology  Demonstrations;  however,  the  actual  sensor  data  was  not  provided 
to  the  Government  and  the  baseline  ordnance  item  details  were  not  released  by 
the  Government  until  after  the  primary  period  of  execution  of  this  investigation. 
Future  data  collected  at  JPG  will  be  provided  to  the  Govenunent  and  baseline 
ordnance  details  will  be  available  following  the  demonstrations.  The  JPG  Phase 
in  baseline  data  have  now  also  been  released.  To  remedy  the  lack  of  readily 
available  multisensor  datsets  to  support  UXO  detection,  discrimination,  and 
identification  research,  ordnance  signature  databases  are  actively  being 
developed  by  Department  of  Defense  agencies  (Army,  Navy,  Defense  Advanced 
Research  Projects  Agency  (DARPA),  and  others),  private  industry,  and 
universities.  The  sections  below  review  the  available  ordnance  signature 
datasets  identified  at  the  initiation  of  this  project  and  efforts  undertaken  as  part 
of  the  project  to  acquire  well-controlled  magnetic  and  acoustic  ordnance 
signatures. 


Development  of  WES  Ordnance 
Signature  Datasets 


Early  in  this  project,  the  requirement  was  established  for  geophysical  signa¬ 
ture  datasets  of  known  ordnance  items  to  validate  phenomenological  modeling 
developments.  The  main  thrust  of  the  work  was  measurement  of  high-accuracy, 
high-resolution  magnetic  signatures  of  inert  ordnance,  using  a  single,  stationary 
magnetometer,  to  support  magnetic  modeling  capability  development.  A 
supplementary  thrust  was  measurement  of  acoustic  signatures  over  ordnance,  to 
investigate  persistent  reports  of  the  existence,  diagnostic  utility,  and  practical 
utility  of  acoustic  resonances  from  ordnance  (e.g.,  Baum  1996).  The  test 
facility  is  located  in  a  wooded  area  on  the  campus  of  Wright  State  University, 
Dayton,  Ohio.  Ordnance  items  in  the  measurement  study  were:  60-mm  mortar 
rounds,  81-mm  mortar  rounds,  90-mm  artillery  shells,  105-mm  artillery  shells 
and  155-mm  artillery  shells. 
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For  the  magnetic  signature  studies  each  ordnance  item  is  placed  in  an  ori- 
entable  holder  made  of  wood  with  non-magnetic  metal  fasteners.  The  center  of 
the  holder  is  orientable  about  a  horizontal  axis  to  allow  positioning  at  any  incli¬ 
nation.  The  holder  is  mounted  on  a  turntable  that  allows  orientation  at  any  azi¬ 
muthal  angle  relative  to  magnetic  north.  The  site  is  a  wooded  area  of  campus;  a 
preliminary  magnetic  survey  determined  the  absence  of  magnetic  anomalies. 
After  the  holder  assembly  is  positioned  and  leveled,  but  before  an  ordnance  item 
is  placed  in  the  holder,  magnetic  readings  are  taken  along  the  track  which  holds 
the  magnetometer.  This  gives  background  profiles  that  show  the  field  to  be 
very  uniform. 

For  the  initial  set  of  magnetic  measurements,  the  magnetometer  sensor  is 
moved  in  10-cm  increments  along  a  track  that  is  oriented  north-south  and  is 
adjusted  to  hold  the  sensor  1.2  m  above  the  center  of  the  ordnance,  except  for 
the  60-mm  mortar  round  for  which  a  smaller  height  is  used.  Measurements  are 
acquired  with  a  proton  precession  magnetometer.  Three  readings  are  averaged 
at  each  position.  A  base  station  is  established  away  from  the  zone  of  influence 
of  the  ordnance  and  readings  are  taken  there  at  least  once  an  hour.  These  read¬ 
ings  are  used  to  remove  the  effects  of  time-variation  of  the  earth's  field  from  the 
measured  data. 

In  the  initial  series  of  TFM  measurements,  the  90-mm  artillery  shell  was 
most  extensively  studied.  For  all  of  the  items,  the  datasets  include  magnetic 
profiles  measured  along  a  profile  across  their  center,  with  measurements  at  10- 
cm  spacing  along  the  profile.  This  was  done  at  the  following  combinations  of 
inclination  angle  (relative  to  horizontal)  and  azimuthal  angle  (in  degrees): 


Inclination 

0 

Azimuth 

0,  45,  and  90 

Inclination 

30 

Azimuth 

0,  45,  and  90 

Inclination 

60 

Azimuth 

0  and  90 

Inclination 

90 

Azimuth 

0 

For  the  90-mm  artillery  shell,  the  dataset  also  includes  these  measurements 
along  north-south  profiles  offset  by  40  cm  and  80  cm  from  the  center  of  the 
shell.  After  base  station  corrections  and  alignment  adjustments,  the  resulting 
curves  show  smooth  variation  of  the  magnetic  field  along  each  profiles.  There 
are  major  differences  in  the  magnitude  and  shape  of  the  profiles  with  ordnance 
orientation,  such  as  illustrated  in  Figure  10  for  measurements  over  a  105-mm 
artillery  projectile.  The  TFM  profiles  are  well-defined  by  the  10-cm  measure¬ 
ments  spacing,  and  the  smooth  profile  curves  indicate  high  accuracy  measure¬ 
ments.  Additional  results  of  the  ordnance  magnetic  signature  measurements  are 
given  in  Appendix  A. 

The  acoustic  studies  are  preliminary  in  nature  and  seek  expression  of  char¬ 
acteristic  acoustic  resonances  that  might  allow  identification  of  various  ordnance 
types.  The  same  ordnance  items  listed  above  are  studied,  with  two  types  of 
measurement  approaches.  In  the  first  type,  a  microphone  is  placed  near  the 
item  and  the  output  of  the  microphone  is  input  to  a  spectrum  analyzer.  The 
ordnance  are  struck  with  a  hard  object  and  the  resulting  sound  spectrum 
recorded.  The  90-mm  and  105-mm  projectiles  produce  characteristic  reso¬ 
nances  at  6.5  kHz  and  5.8  kHz,  respectively.  However,  the  other  three  sizes 
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All  data  corrected  to  common  base  station  (B.S.)  measurement:  105-mm  Projectile 
0-cm  offset  profiles;  all  orientations;  background  measurements  subtracted  (subtr) 


Chapter  3  Data  Sources  for  the  Investigation 


23 


Figure  10.  Example  from  WES  ordnance  signatures  datasets;  TFM  measurements  over  center  of  105-mm  projectile,  1 .2  m  below  sensor, 
and  for  nine  orientations 


do  not  produce  obvious  resonances  that  are  strong  and  consistent  between  sam¬ 
ples  and  relatively  independent  of  the  location  of  the  strike. 

For  the  second  type  of  measurement,  an  item  is  placed  in  a  home-made 
acoustic  enclosure.  A  microphone  is  attached  to  the  backside  of  the  item  and  a 
speaker  one  meter  away  is  swept  through  a  range  of  frequencies  from  3  kHz  to 
10  kHz.  Again,  the  strongest,  most  consistent  resonances  are  at  the  same  fre¬ 
quencies  as  before  for  the  90  mm  and  105  mm  projectiles.  Subsequently,  the 
items  were  covered  to  various  depths  in  wet  and  dry  sand;  with  the  perhaps 
intuitive  result  that  the  strength  of  the  resonances  progressively  decrease  as 
more  of  the  item  was  covered. 

Impulsive  measurements  in  a  container  filled  with  sand  or  water  likewise 
indicate  that  the  resonance  strengths  progressively  decrease  as  the  items  are 
more  deeply  buried.  Preliminary  and  tentative  conclusions  of  the  acoustic  tests, 
for  the  limited  suite  of  ordnance  tested,  are  that  two  ordnance  types  have 
characteristic  resonances,  but  the  resonances  are  too  damped  on  burial  to  pro¬ 
vide  a  viable  identification  technique  for  buried  ordnance.  A  subsequent  report 
will  present  the  acoustic  results  in  detail  as  well  as  all  measured  magnetic 
signatures. 


Naval  Research  Laboratory  (NRL) 
Ordnance  Signature  Library 


The  Multisensor  Towed  Array  Detection  System  (MTADS)  is  the  culmina¬ 
tion  of  at  least  a  decade  of  engineering  development  (McDonald  and  Robertson 
1996).  The  MTADS  has  overcome  many  of  the  difficulties  of  multisensor  inte¬ 
gration  on  a  towed  platform.  MTADS  utilizes  separate  towed  platforms  to 
acquire  high-resolution  passive  magnetic  (eight  TFM  sensors  or  four  TFG 
sensor  sets)  and  active  electromagnetic  (three  overlapping  TDEM  sensors)  mea¬ 
surements.  A  differential  global  positioning  system  (DGPS)  is  used  for  primary 
navigation  in  the  field.  To  support  modeling  and  detection  and  discrimination 
algorithm  development,  an  ordnance  signature  library  is  being  developed  using 
MTADS  (Nelson,  McDonald,  and  Robertson  1997,  Barrow  et  al.  1997). 

For  development  of  the  NRL  ordnance  signature  library,  MTADS  sensors 
are  used  to  acquire  TFM,  TFG  and  TDEM  (Geonics  EM61)  datasets  over  inert 
ordnance  placed  on  the  surface,  in  a  1  m  deep  hole,  and  in  a  7  m  deep  well. 
Special  jigs  hold  the  ordnance  items  at  predetermined  depths  and  orientations. 
Tables  3,  4  and  5  summarize  the  datasets  from  the  signature  library  that  were 
available  from  NRL  for  use  in  the  present  investigation.  As  an  example,  TFM 
and  TDEM  datasets  from  the  NRL  signature  library  for  a  105  mm  projectile 
(depth  =  0.49  m;  inclination  =  0;  azimuth  =  90)  are  shown  in  Figure  11.  The 
datasets  in  Figure  1 1  are  shown  as  images,  i.e.,  color-filled  contour  maps.  The 
NRL  ordnance  library  is  used  for  phenomenology  modeling  validation 
(primarily  TFM  modeling)  and  for  parameter  cross-plotting  investigations. 
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Table  3 

Ordnance  Items,  Depths  and  Orientations,  for  the  TFM  Data  Sets, 

NRL  Ordnance  Signatures  Library 

Item 

Depths 

Azimuth 

Inclination 

20  mm  projectile 

Surface 

o 

o 

(D 

O 

o 

0® 

30  mm  projectile 

Surface 

0 

o 

O) 

0 

o 

0® 

M42  grenade 

Surface,  1 5  cm 

0%  90** 

0® 

M46  submunition 

Surface,  1 5  cm 

0°,  90® 

0® 

60  mm  mortar 

0.25,  0.5  m 

45®  steps 

45®  steps 

81  mm  mortar 

0.5,  0.75,  1  m 

45®  steps 

45®  steps 

105  mm  projectile 

0.5,  0.75,  1  m 

45®  steps 

45®  steps 

5-in.  rocket 

1,  1.5  m 

45®  steps 

45®  steps 

Table  4 

Ordnance  Items,  Depths  and  Orientations,  for  the  TFG  and  EM61  Data  Sets, 
NRL  Ordnance  Signatures  Library 

Item 

Gradiometer  Survey 

EM-61  Survey  || 

Depths 

Azimuth 

Inclination 

Depths 

Azimuth  and 
Inclination 

20  mm  projectile 

Surface 

o 

o 

<D 

O 

o 

0® 

Surface 

0®,  90® 

30  mm  projectile 

Surface 

0 

o 

0) 

o 

o 

0® 

Surface 

0®,  90® 

M42  grenade 

Surface,  15  cm 

0®,  90® 

0® 

Surface,  1 5  cm 

o 

o 

0> 

o 

o 

M46  submunition 

Surface,  1 5  cm 

o 

o 

(O 

o 

0 

0® 

Surface,  15  cm 

o 

0 

(O 

o 

0 

60  mm  mortar 

0.25  m 

45®  steps 

45®  steps 

0.25,  0.5,  0.75,  1  m 

90®  steps 

81  mm  mortar 

0.5,  0.75  m 

45®  steps 

45®  steps 

0.5,  0.75,  1  m 

90®  steps 

1 05  mm  projectile 

0.5,  0.75  m 

45®  steps 

45®  steps 

0.5,  0.75,1,  1.25  m 

90°  steps 

5-in.  rocket 

1  m 

45®  steps 

45®  steps 

0.5,  1,  1.5  m 

90°  steps 

155  mm  projectile 

1.5,  2  m 

90®  steps 

DARPA  Background  Clutter  Program 


In  addition  to  the  preceding  ordnance  signature  libraries  developed  in  well- 
controlled,  “laboratory-type”  settings,  two  additional  sources  of  well- 
controlled,  multisensor  datasets  acquired  in  settings  more  “real  world”  are  uti¬ 
lized  in  the  investigation.  The  DARPA-sponsored  Background  Clutter  Data 
Collection  Experiment  was  designed  to  address  a  void  in  multisensor  data  col¬ 
lected  specifically  to  characterize  the  natural  background  and  manmade  clutter 
at  sites  where  buried  UXO  and  landmines  are  typically  found  (George  and 
Altshuler  1997).  Within  the  limited  scope  of  the  effort,  a  variety  of  geologic 
settings  was  accommodated  at  four  test  sites  at  two  geographic  locations:  Fort 
A.  P.  Hill,  Virginia— two  “wet,  sand”  sites.  Firing  Point  20  (FP20)and  Firing 
Point  22  (FP22);  Fort  Carson,  Colorado— a  “dry,  sand”  site,  Turkey  Creek 


Chapter  3  Data  Sources  for  the  Investigation 


25 


Table  5 

Mass  and  Dimensions  of  Ordnance  Items,  NRL  Ordnance  Signatures 
Library 

Item 

Length  (cm) 

Diameter  (cm) 

Mass  (kg) 

60-mm  mortar 

35.6 

6.0 

1.0 

81 -mm  mortar 

42 

8.1 

2.3 

2.75-in.  rocket 

33 

7.0 

4.1 

4.2-in.  mortar 

45 

10.7 

8.6 

105-mm  projectile 

39 

10.5 

9.5 

5-in.  rocket 

49 

12.7 

20 

1 55-mm  projectile 

62 

15.5 

25.4 

Mk81  bomb  (250  lb) 

23.0 

56.7 

8-in.  projectile 

78.7 

20.3 

73.9 

Mk82  bomb  (500  lb) 

156 

28.0 

113 

Mk83  bomb  (1000  lb) 

191.5 

36.6 

227 

Site,  and  a  “dry,  clay”  site,  Seabee  Site.  Each  test  site  was  a  standard  125  m  x 
100  m  with  a  central  lOO-m  x  100-m  backgrounds  area  (Center  Square).  The 
only  emplaced  objects  in  the  Center  Square  were  five  sets  of  registration  targets, 
each  set  consisting  of  a  buried  iron  sphere  and  an  aluminum  plate  spaced  4  m 
apart  (an  additional  aluminum  plate  was  placed  on  the  surface  midway  between 
the  buried  sphere  and  plate  for  infrared  ^R)  data  collection).  On  two  sides  of 
the  sites,  target  sidebars  were  established,  containing  a  variety  of  buried  inert 
ordnance,  inert  landmines,  and  other  objects.  Inert  ordnance  targets  consist  of 
60  mm  and  81  mm  mortars  and  105  mm  and  152  mm  projectiles.  The  general 
site  layout  is  shown  in  Figure  12.  Details  regarding  target  locations  in  the 
sidebars  has  not  been  released,  except  for  the  region  noted  as  the  Calibration 
Area  in  Figure  12  (shaded).  Individual  target  locations  and  descriptions  for  the 
Calibration  Area  are  shown  in  Figure  13  (George  and  Altshuler  1997). 

Seven  contractors  obtained  high-resolution  multisensor  data  over  the  test 
sites.  Although  data  density  over  the  sites  varied  from  one  sensor  system  to 
another,  typical  data  density  is  10  to  20  cm  along  track,  with  25-  to  50-cm  track 
spacing.  Data  acquired  by  Geophex,  Geometries,  and  Parsons  were  provided 
directly  to  WES  and  also  to  DARPA.  All  data  acquired  under  the  Backgrounds 
Data  Collection  Experiment  are  available  on  approximately  35  CD-ROM’s,  that 
are  distributed  for  DARPA  by  Walcoff  and  Associates  (George  and  Altshuler 
1997).  Although  direct  correlation  of  processed  data  with  the  baseline  target 
cannot  be  presented  in  this  report  (except  for  the  Calibration  Area),  the  datasets 
are  useful  for  illustrating  target  detection  in  general  and  for  illustrations  of  data 
management  and  of  “image  processing”  to  enhance  anomalies.  The  following 
tabulation  briefly  summarizes  the  types  of  sensor  data  collected  at  the  test  sites: 
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Sensor/Method 

Contractor 

Details 

Magnetometry 

Geometries 

Six  Geometries  G-858  cesium  vapor  magnetometers. 
Mounted  on  cart;  configured  to  give  high-density 
total  field  magnetic  and  horizontal  and  vertical  mag¬ 
netic  gradient. 

Geophex 

Two-element  Geometries  G-858  total  field  magne¬ 
tometer  system. 

EMI 

Coleman  Research 
Corp 

Three-element  Geonics  EM61  (TDEM)  array  (1-m 

Tx/Rx) 

Parsons  Engineering 
Science,  Inc, 

Custom  built,  0.5  m  Tx  and  Rx,  Geonics  EM61 
(TDEM). 

Prototype  multichannel,  multicomponent  EM61-3D, 
with  1  m  Tx  and  3-orthoqonal  axis  0.5  m  Rx  (TDEM) 

Science  Applications 
International  Corp. 

Array  of  eight  overlapping  Schiebel  EMI  coils  (TDEM). 

Geophex 

Prototype  Multifrequency  FDEM  system,  operated  at 
4,050  Hz  and  12,270  Hz. 

GPR 

Coleman  Research 
Corp. 

GPR  array  (2  Tx,  3  Rx).  Frequency  domain;  2  MHZ 
stepped-frequency  over  range  100-610  MHZ.  Addi¬ 
tional  Tx-Rx  with  10  MHZ  stepped-frequency  over 
range  1,000  -  4,000  MHZ.  Vehicle  towed. 

— - - — 1 

Geocenters,  Inc. 

"Focused  Array  Radar"  (4  Tx,  4  Rx).  Frequency 
domain,  swept  frequency  over  range  700- 
1.300  MHZ. 

Lawrence  Livermore 
National 

Laboratories 

Prototype  GPR.  Frequency  domain  over  range  5  MHZ 
-  18  GHz.  with  up  to  801  steps.  Operated  in 
monostatic,  cross-polarized  mode. 

IR 

Geo-Centers,  Inc. 

Forward-looking  system. 

Science  Applications 
International,  Inc. 

Down-looking  system. 

Datasets  from  MTADS  Surveys  at  the  Magnetic 
Test  Range,  Twenty-nine  Palms,  California 


Another  well-controlled,  “real  world”  multisensor  dataset  source  is  Ae 
MTADS  survey  results  at  the  Magnetic  Test  Range  (MTR)  at  Twenty-nine 
Palms,  California  (McDonald  et  al.  1997).  The  MTR  is  approximately  220  m 
X  150  m  (about  8  acres)  in  size.  Table  6  lists  the  ordnance,  number  of  items, 
and  associated  burial  depth  range  of  the  70  ordnance  items  buried  at  the  site 
(from  McDonald  et  al.  1997).  Three  complete  surveys  of  the  site  are  available: 
TFM,  TFG,  TDEM.  For  the  TFM  survey,  the  sensors  were  set  0.25  m  above 
the  surface,  and  data  are  acquired  at  approximately  0.06  m  along  track,  with 
0.25-m  track  spacing;  while  for  the  TFG  survey,  the  sensors  were  reconfigured 
with  sensor  pairs  at  0.40  m  and  0.95  m  above  the  surface  and  0.5-m  track 
spacing.  The  lower  transmit/receive  coils  of  the  TDEM  array  (3  EMbl's)  were 
0.40  m  above  the  surface,  and  data  are  acquired  approximately  0. 15  m  along 
track,  with  0.5-ni  track  spacing. 
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Figure  1 1 .  NRL  ordnance  signature  library  example;  TFM  and  IDEM  data  images  for  105-mm  projectile,  depth  =  0.49  m,  inclination  =  0  deg, 
azimuth  =  90  deg 


Registration 
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Figure  12.  Basic  layout  of  the  DARPA  1 -hectare  sites 


3' 


Figure  13.  Calibration  area  layout  for  the  DARPA  1 -hectare  sites 


Table  6 

Ordnance  Inventory  at  the  Twenty-nine  Palms  Mi 
(MTR) 

agnetic  Test  Range 

Ordnance 

Number  of  Items 

Range  of  Depths  (m) 

60-mm  Mortar 

10 

0.15-0.46 

81 -mm  Mortar 

7 

0.46-0.76 

105-mm  Projectile 

10 

0.46-1.10 

155*mm  Projectile 

10 

0.61-1.22 

8-in,  Projectile 

10 

1.83-2.74 

Mk  81  Bomb 

10 

1.43-3.11 

Mk  82  Bomb 

10 

1 .22-4.42 

Mk 117  Bomb 

1 

3.96 

Mk  83  Bomb 

1 

5.09 

Mk  84  Bomb 

1 

4.88 
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Data  Management  and 
Analysis  Methods  for 
Detection  Enhancement 


The  importance  of  data  management  and  analysis  tools  is  emphasized  in 
Chapter  2,  and  selected  examples  are  presented  here  using  data  from  the 
DARPA  sites.  As  mentioned  previously,  there  are  three  approaches  to  meeting 
the  data  management  and  analysis  requirements  for  multiple  (multisensor)  data¬ 
sets;  (1)  develop  special  purpose  data  acquisition,  management  and  analysis 
software  for  a  specific  system  (e.g.,  McDonald  and  Robertson  1996);  (2)  utilize 
commercially  available  general  purpose  graphics  software  that  has  some  GIS- 
type  functionality,  such  as  simultaneous  display  of  two  or  more  maps  and  data 
linkages  between  the  maps  (e.g.,  Geosoft  1996);  (3)  use  a  full  functioned  GIS 
for  data  display,  mapping,  relational  links,  and  analyses. 


GIS  Data  Management  Examples 


Data  for  the  DARPA  Fort  Carson  sites,  including  site  characterization 
results  and  contractor  sensor  datasets,  were  input  to  the  ARCVIEW  *  GIS  (ESRI 
1997).  The  DARPA  site  characterization  data,  acquired  prior  to  any  burial 
activities  at  the  sites,  included  detailed  mapping  of  topography,  vegetation,  and 
surface  features  (e.g.,  rocks,  animal  burrows,  metallic  objects,  fences,  etc.),  as 
well  as  geophysical  surveys,  soil  sampling,  and  soil  classifications  and 
laboratory  EM  property  measurements.  The  first  example.  Figure  14,  is  the 
topographic  map  overlying  the  key  site  features  map  for  the  Turkey  Creek  Site, 
Fort  Carson.  Clicking  on  a  site  feature  brings  up  an  information  box  describing 
all  details  of  the  feature.  Tbe  “x’s*  indicate  surface  site  features,  e.g.,  rocks  or 
cultural  features,  and  the  solid  circles  in  Figure  14  are  the  locations  of  the 
releasable  buried  targets.  The  display  in  Figure  14  is  an  important  benchmark 
map  for  correlation  with  all  site  characterization  geophysical  surveys  and  with 
contractor  geophysical  sensor  datasets.  Topographic  features,  surface 
vegetation  and  roots,  rocks,  and  metallic  cultural  features  can  all  produce 
geophysical  anomalies  which  may  contribute  to  the  false  alarm  rate,  affecting 
UXO  detection  and  discrimination.  Another  example  illustrates  a 
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Figure  14. 


simple  query  of  two  site  characterization  maps  (  TFM  and  EM31^)  for  the  Fir¬ 
ing  Point  20  Site,  Fort  A.  P.  Hill:  The  query  computes  an  average  value  for  the 
geophysical  measurements  defining  the  TFM  map  and  the  EM31  (FDEM)  map. 
The  query  then  produces  a  map.  Figure  15,  which  indicates  all  areas  where  both 
the  TFM  and  EM31  measurements  are  greater  than  the  average  values.  Two 
linear  features  are  defined  in  Figure  15,  trending  from  southwest  to  northeast  at 
the  top  of  the  map.  These  linear  features  are  prominent  on  many  of  the 
subsequent  TFM,  TDEM,  and  FDEM  datasets  acquired  at  the  site.  These 
features  were  subsequently  verified  as  caused  by  buried  communications  cables. 
In  general  >the  coincident  anomaly  areas  defined  in  Figure  15  indicate  locations 
where  sensor  datasets  may  be  more  complicated  to  interpret  and  are  possible 
sources  of  false  alarms. 

The  GIS  provides  a  versatile  environment  for  display  of  sensor  “images.” 

For  example.  Figure  16  is  an  image  or  map  of  TDEM  measurements  (EM61, 

0.5  m  Tx-Rx)  for  the  Turkey  Creek  Site  at  Fort  Carson.  Color  values  in  the 
map  are  based  on  multiples  of  the  standard  deviation  above  or  below  the  mean 
value  for  the  dataset.  The  map  indicates  localized,  high-resolution  anomalies 
indicative  of  known  buried  metallic  objects  and  some  localized  anomalies  of 
unknown  origin  (false  alarms).  All  five  of  the  registration  target  sets  are  clearly 
indicated.  It  is  a  simple  procedure  to  enlarge  selected  areas  of  the  map  for  more 
detailed  examination  of  anomalies.  The  enlarged  example  in  Figure  16  includes 
two  sets  of  the  registration  targets,  with  the  anomalies  identified  with  the  actual 
buried  objects  for  one  of  the  sets.  The  aluminum  plate  and  iron  sphere  of  the 
registration  target  sets  are  spatially  separated  by  4  m  (center  to  center);  the 
anomalies  caused  by  the  targets  are  well  resolved.  There  are  several  small 
magnitude  anomalies  that  could  be  due  to  unknown  buried  metallic  objects  (note 
the  small  spatial  extent  anomalies  just  to  the  right  of  the  lower  right  registration 
targets  in  the  enlarged  view). 

Another  useful  tool  with  the  GIS  is  the  capability  of  easily  performing  data 
transformations  on  one  or  more  map  layers.  The  background  noise  and  clutter 
in  Figure  16  is  evident  and  was  discussed  in  the  preceding  paragraph.  A  useful 
technique  for  suppressing  background  and  clutter  is  setting  a  display  threshold. 
Figure  17  is  a  map  of  all  values  in  Figure  16  that  exceed  the  mean  plus  one 
standard  deviation  (displayed  in  red).  In  the  central  area,  the  background  and 
clutter  contributing  to  potential  false  alarms  is  virtually  eliminated. 


Image  Processing  of  Multisensor 
Data  for  Detection  Enhancement 


A  conunercially  available  graphics  software  package  that  has  some  GIS 
functionality,  including  simultaneous  display  of  multiple  co-registered  data 
maps  and  some  relational  links  between  maps,  is  OASIS  montaj  (Geosoft 


‘  The  EM31  is  a  single  frequency  (9.6  kHz;  Tx-Rx  spacing  =  3.7  m)  EMI  system  manufactured 
by  Geonics,  Inc.,  Mississauga,  Ontario,  Canada. 
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Figure  1 5.  Results  of  GIS  query  of  WES  site  characterization  data  for  the  Firing  Point  22  Site,  Fort  A.P.  Hill.  The  query  computes  average  values  for  the 
TFM  and  EM3 1  (conductivity)  maps  and  produces  a  derived  map  showing  all  areas  where  both  TFM  and  EM3 1  values  are  above  the  averages 
(in  black).  Vertical  black  bars  or  strips  are  areas  of  no  data.  Distance  in  meters 
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Figure  16.  GIS  map  of  EM61  (0.5«m  Tx-Rx)  for  the  Turkey  Creek  Site,  Fort  Carson.  Zoom  capability  shows  close-up  of  anomalies  caused  by  two 
of  the  sets  of  registration  targets  in  the  Center  Square.  Vertical  black  bars  indicate  no  data  acquired.  Grid  line  spacing  —  5  m  x  5  m 
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deviation.  Distance  is  in  meters;  north  is  to  top  of  figure 


1996).  Figure  18  is  an  example  of  a  side-by-side  presentation  of  TFM  analytic 
signal  (see  Chapter  5  for  a  definition  of  the  magnetic  analytic  signal)  and  TDEM 
(EM61)  image  maps  for  the  Firing  Point  20  Site,  Fort  A.  P.  Hill.  The  images 
are  a  shade  of  gray  version  of  maps  originally  in  color;  a  single  and  different 
color  for  each  image,  with  color  intensity  increasing  as  data  magnitude 
increases.  There  is  considerable  clutter  noise  in  both  images;  although  some 
small  magnetic  clutter  signals  are  suppressed  in  the  analytic  signal  repre¬ 
sentation.  Both  of  these  maps  can  be  treated  as  images,  with  pixel  size  defined 
by  the  basic  gridding  and  data  acquisition  spacings.  The  original  color  images 
were  combined  to  form  a  color-merged  composite  or  false  color  image.  Fig¬ 
ure  19  is  the  shade  of  gray  version  of  the  color  composite  image,  where  increas¬ 
ing  brightness  (i.e.,  black  to  white)  indicates  increasing  intensity  of  both 
original  color  images.  For  example,  a  localized  “w^ite  anomaly”  in  the  com¬ 
posite  image  is  a  target  or  feature  with  both  high  magnetic  analytic  signal  mag¬ 
nitude  and  high  EM61  magnitude.  Many  of  the  potential  false  alarm  targets  in 
the  two  original  images  are  suppressed,  many  localized  anomalies  representing 
potential  targets  of  interest  are  evident  and  enhanced,  and  several  linear  clutter 
features  are  well-defined. 
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Figure  18.  Shade  of  gray  versions  of  color  image  maps  of  magnetic  analytic  signal  and  IDEM  (EM61)  for  Firing  Point  20,  Fort  A.P.  Hill,  Virginia 
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Figure  1 9.  Shade  of  gray  version  of  color-merged  composite  image  of  magnetic  analytic  signal  and  IDEM  (EM61)  images  for  Firing  Point  20, 
FortA.P.  Hill 


5  Synopsis  and  Examples  of 
Current  UXO  Detection 
Analyses  Capability 


Introduction 


While  there  is  predominantly  only  rudimentary,  empirical  multisensor 
(multimethod)  integration  in  fielded  UXO  detection  systems,  detection  and 
depth  interpretation  capability  from  single-method  datasets  is  actually  quite 
good,  particularly  for  fielded  TFM  and  TDEM  systems.  The  primary 
problem,  as  discussed  earlier,  is  that  there  is  not  a  significant  discrimination 
capability  based  on  single  method  data.  Although  there  is  some  automated 
anomaly  selection  capability,  many  of  the  interpretation  approaches  require 
manual  (interpreter)  anomaly  picking  or  selection  (often  consisting  of  “boxing- 
in”  the  anomaly).  Interpretation  of  a  selected  TFM  or  TDEM  anomaly 
proceeds  by  applying  an  empirical  or  analytical  model  to  the  anomaly  character¬ 
istics.  Currently,  model  based  interpretation  generally  relies  on  simple 
geometric  models,  such  as  spheres  or  cylinders.  For  example,  most  available 
TFM  interpretation  approaches  use  a  simple  dipole  model  (equivalent  to  a  uni¬ 
formly  magnetized  sphere)  to  fit  anomaly  data.  From  the  simple  dipole  model 
TFM  anomaly  interpretation,  the  position  and  depth  coordinates  and  dipole 
magnetic  moment  (magnitude  and  orientation)  of  a  best-fitting  magnetic  dipole 
are  determined.  The  best-fitting  simple  dipole  magnetic  moment  magnitude 
gives  an  equivalent  sphere  size  (radius  and  ferrous  mass;  e.g.,  Pennella  1982), 
which  correlates  to  ordnance  diameter  (Barrow  et  al.  1997). 

Since  the  field  of  a  simple  magnetic  dipole  is  equivalent  to  the  field  of  a 
uniformly  magnetized  sphere,  the  earth’s  field-induced  magnetization  for  this 
model  (in  the  absence  of  remnant  magnetization)  will  be  in  the  direction  of  the 
earth’s  field.  However,  the  simple  dipole  model  fit  to  measured  data  over 
UXO  are  frequently  not  in  the  direction  of  the  earth’s  field,  because  ordnance 
items  are  not  uniformly  magnetized  spheres.  The  effects  of  the  elongated 
shape  of  UXO  on  induced  magnetic  moments  in  the  earth’s  field  are  studied  by 
McFee  and  Das  (1990)  and  Altshuler  (1996)  using  a  prolate  spheroid  model. 
Barrow  et  al.  (1997)  investigate  the  ambiguities  resulting  from  variations  in 
best-fitting  simple  dipoles  due  solely  to  changes  in  orientation  of  the  UXO. 
Barrow  et  al.  also  investigate  the  capability  of  the  prolate  spheroid  model  to 
replicate  the  variation  in  best-fitting  dipole  (magnitude  and  orientation)  as  a 
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function  of  size  and  orientation  of  the  UXO.  The  prolate  spheroid  model 
predictions  of  the  variation  differ  only  in  small  details,  likely  resulting  from 
slight  symmetry  breaking  along  the  long  axis  of  real  ordnance.  Application  of 
the  simple  magnetic  dipole  model  to  interpretation  of  TFM  anomalies  achieves 
good  performance  for  determining  position  and  depth,  with  typical  errors  of 
±  10  to  20  percent.  It  cannot  be  used  to  determine  the  size  (length  and 
diameter)  and  orientation  (Barrow  et  al.  1997)  uniquely  or  accurately. 

Similarly  for  the  TDEM  analysis,  a  spherical  model  is  often  used  to  estimate 
depth  and  size  (Barrow  et  al.  1997).  For  the  EM61  TDEM  system,  a  simple 
algebraic  formula  allows  depth  estimates  based  on  the  ratio  of  an  upper  to  a 
lower  receiver  coil  signals  ^ollyer,  Racic,  and  Butler  1997).  The  simple 
dipole  model  for  magnetic  interpretation  and  the  sphere  model  for  TDEM 
interpretation  cannot  account  for  effects  of  ordnance  elongation  on  the  induced 
anomalies.  The  breakdown  of  the  simple  dipole  or  sphere  models  for 
interpreting  anomalies  becomes  most  pronounced  for  objects  buried  less  than 
(1)  the  ordnance  length  in  depth  for  TFM  or  (2)  one  to  two  times  the  effective 
transmitter  radius  for  TDEM.  A  great  advantage  to  the  simple  model 
interpretations  is  computational  speed;  the  interpretations  can  be  accomplished 
in  near-real  time.  Two  examples  that  illustrate  current  capabilities  and  a 
prototype  multisensor  integration  algorithm  are  presented  below. 


Fits  to  Twentynine 
Palms  MTADS  Datasets 


As  discussed  previously,  discrimination  capability  for  reduction  of  false 
alarms  is  key  to  reducing  UXO  cleanup  time  and  cost.  Accurate  burial  depth 
estimates  for  UXO  are  also  extremely  important  for  UXO  remediation/cleanup 
time  and  cost  estimates  and  for  safety  considerations  during  excavation  of  UXO 
items.  Based  on  considerable  experience  at  test  sites  and  at  live  sites,  i.e., 
predictions  based  on  interpretations  of  geophysical  anomalies  compared  to 
actual  excavation  details  of  the  features  producing  the  anomalies,  there  is  an 
understanding  of  the  accuracy  and  reliability  of  simple  model  predictions.  Sim¬ 
ple  magnetic  dipole  model  fits  of  ordnance  depth  and  size  to  TFM  data  are  typi¬ 
cally  good  to  ±  10  to  15  percent  for  location,  ±  20  percent  for  depth,  and 
±  25  percent  for  size  (e.g..  Bell  1997). 

The  results  of  the  MTADS  surveys  at  Twentynine  Palms  successfully  demon¬ 
strated  detection  capability  in  a  real  world  setting,  with  probabilities  of 
detection  (for  a  joint  or  “fused”  analysis)  ranging  from  73  percent  for  0.5-m 
critical  radius  to  94  percent  for  2.0-m  critical  radius,  where  a  detection  occurs 
when  a  declared  ordnance  location  is  within  the  specified  critical  radius  of  a 
true  ordnance  location.  The  false  alarm  ratio  metric  (number  of  false  alarms 
per  ordnance  item  detected)  ranged  from  1.9  for  individual  analyses  of  magnetic 
or  TDEM  anomalies  to  3.0  for  a  joint  analysis.  Joint  analysis  of  the  three 
MTADS  datasets  consisted  of  visual  correlation,  with  selection  and  deselection 
based  on  anomaly  coincidence,  type,  and  predicted  size  and  depth 
considerations.  The  joint  analysis  resulted  in  three  additional  valid  ordnance 
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target  detections  (for  a  2-m  critical  radius),  raising  the  probability  of  detection 
to  94  percent  but  significantly  increasing  the  number  of  false  alarms.  The  joint 
analysis  results  here  are  consistent  with  the  experience  at  JPG.  Joint  analysis  as 
currently  practiced  desirably  increases  the  probabilities  of  detection  but  also 
undesirably  increases  the  number  of  false  alarms. 

Due  to  the  wide  and  realistic  depth  range  of  ordnance  at  the  Twentynine 
Palms  site,  interpreted  or  fitted  depths  provide  a  good  assessment  of  depth 
estimation  capabilities  for  simple  model  fits.  Comparisons  of  fitted  depth  to 
true  depth  for  some  of  the  ordnance  items  detected  by  each  of  the  three  MTADS 
sensor  suites  are  shown  in  Figure  20.  The  fitted  depths  in  Figure  20  generally 
cluster  about  the  dashed,  “perfect-fit”  line.  There  are  a  few  outliers  in  the 
gradiometer  fitted  depths,  and  the  TDEM  (EM-61)  fitted  depths  deviate 
considerably  for  true  depths  comparable  to  and  less  than  the  transmitter  equiva¬ 
lent  radius  (approximately  0.6  m).  Figure  21  contains  a  comparison  of  fitted 
size  from  the  TDEM  data  versus  true  caliber  (diameter)  and  against  an  equiva¬ 
lent  sphere  radius.  For  a  given  ordnance  size  (caliber  or  equivalent  sphere 
radius),  the  mean  of  the  fit  sizes  correlates  well  with  the  true  size.  However, 
the  spread  of  TDEM  fit  sizes  varies  by  factors  of  typically  2  to  3  for  a  given 
true  ordnance  size.  The  spread  in  computed  depths  and  sizes  can  be  partially 
explained  by  the  fact  that  the  simple  spherical  model  cannot  account  for  varia¬ 
tions  in  ordnance  orientation,  which  can  have  significant  effect  on  TDEM 
anomaly  signature. 


Simple  Model  Fits  to  Magnetic 
Data  from  DARPA  Site 


A  small  window  (5  m  x  50  m)  of  TFM  data  from  the  DARPA  Seabee  Site, 
Fort  Carson,  Colorado,  is  shown  in  Figure  22.  The  Fort  Carson  sites  are 
“quiet”  magnetically,  so  that  anomalies  from  shallow,  buried,  ferrous  objects 
are  pronounced.  In  Figure  22,  five  prominent  TFM  anomalies  are  indicated 
(enclosed)  by  boxes.  A  commercially  available  inversion  program,  MagAID, 
developed  by  AETC,  Inc.,  and  marketed  by  Geometries,  Inc.  (AETC- 
Geometrics  1996),  fits  simple  dipole  models  to  measured  TFM  anomaly  signa¬ 
tures.  MagAID  determines  the  set  of  simple  magnetic  dipole  parameters  which 
best  fit  the  measured  data  by  minimizing  the  Chi-squared  statistic  between  the 
measured  data  and  model  predictions.  The  program  uses  a  modified  gradient 
search  technique  through  the  parameter  space  (Marquardt  1963;  DeProspo  and 
DiMarco  1996).  Applying  the  MagAID  analysis  procedure  to  the  five  selected 
anomalies  in  Figure  22  results  in  the  interpreted  parameters  in  Table  7. 

Twenty  TFM  anomalies  caused  by  buried  objects  (some  known  objects  and 
some  unknown  objects)  at  various  locations  on  the  Seabee  Site  were  inverted  to 
obtain  simple  dipole  fit  model  parameters.  While  the  actual  known  object 
identities,  locations,  and  depths  cannot  be  released,  the  average  errors  of  the 
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MTADS  Depth  Fits  at  Twentynine  Palms 


California  (McDonald  et  all997) 


Figure  21 .  Correlation  of  fitted  size  from  analyses  of  the  Twentynine  Palms  IDEM  data  as  a  function  of 
(a)  caliber  (top)  and  (b)  equivalent  sphere  radius  (bottom) 
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Seabee  Site 
Fort  Carson,  Colorado 
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Figure  22.  Selected  window  (5  m  x  50  m)  of  TFM  data,  Seabee  Site,  Fort  Carson.  Five  prominent  anomalies  selected  and  indicated  by  the  boxed  areas 


Table  7 

Interpretation  of  TFIV 

Anomalies,  Seabee  Site, 

Fort  Carson 

Object 

No. 

Interpreted 

Object  Location 
-  X,  m 

Interpreted 

Object  Location 
-  Y.  m 

Interpreted 

Object 

Depth,  m 

Interpreted 
Object  Size, 
mm’ 

Rt 

Quality 

U^.0) 

6 

122.1 

53.8 

0.60 

112 

0.77 

7 

122.7 

73.9 

0.77 

127 

0.92 

8 

122.1 

80.0 

0.43 

80 

0.99 

9 

122.8 

87.6 

0.83 

110 

0.84 

10 

123.2 

96.4 

1.29 

120 

0.88 

'  TheO 
anomaly 
equivalei 

bject  Size  is  the  probable  ordnance  diameter,  assuming  the  object  causing  the 
is  ordnance.  Optionally  the  Object  Size  can  be  output  as  the  radius  of  a  dipole 
nt  sphere. 

twenty  interpretations  that  correspond  to  known  objects  gives  insight  to  the 
inversion  capabilities  of  the  program: 

Average  Location  Error  -  0.3  m 

Average  Depth  Error  (percent  of  actual)  -  45  percent 

Average  Size  Error  (percent  of  actual)  -  14  percent 

The  average  location  error  is  approximately  the  same  as  the  average  measure¬ 
ment  spacing;  the  measurement  grid  is  approximately  10  to  20  cm  along  lines 
and  25  to  50  cm  between  lines.  The  large  average  depth  error  of  the 
interpretations,  as  a  percentage  of  actual  depth,  is  due  to  the  fact  that  actual 
depths  of  buried  objects  at  the  DARPA  sites  are  small.  There  is  a  strong 
correlation  of  computed  depths  with  actual  depths. 

Conclusions  from  the  error  results  for  the  MTADS  Twentynine  Palms  sur¬ 
veys  and  the  error  results  listed  above  for  the  DARPA  TFM  target  fits  are  con¬ 
sistent  with  the  extensive  results  from  the  recently  reported  MTADS  surveys  at 
the  Badlands  Bombing  Range  (McDonald  et  al.  1997).  A  large  number  of  tar¬ 
gets  were  dug  up  after  TFM  and/or  TDEM  anomaly  detection,  selection,  and 
interpretation.  The  average  error  of  target  locations  (interpreted  versus  dug 
target  location)  was  12  cm,  and  95  percent  of  all  targets  were  located  within 
29  cm.  The  predicted  versus  actual  target  depths  for  all  dug  targets  are  strongly 
correlated  but  show  considerable  variation,  particularly  for  shallow  targets. 
Similar  strong  correlation  between  predicted  size  and  actual  size  exists,  but 
considerable  variation  exists  also.  As  in  previous  examples,  multisensor  data 
integration  for  the  MTADS  Badlands  Bombing  Range  surveys  was  accomplished 
by  visually  correlating  TFM  and  TDEM  anomalies. 
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Proposed  Multisensor 
Data  Integration  Algorithm 


The  examples  of  multisensor  data  integration  discussed  to  this  point  are 
predominantly  empirical,  involving  side-by-side  visualization  of  datasets  and 
selection  of  spatially  coincident  anomalies  on  two  or  more  datasets,  e.g.,  the 
MTADS  examples  and  the  example  in  Figure  18.  An  intermediate  stage  of 
multisensor  data  integration,  between  the  visual  correlation  procedure  and  true 
joint  inversion  of  multiple  datasets,  is  to  invert  one  dataset  using  selected  results 
from  the  inversion  of,  or  parameter  extraction  from,  a  second  dataset;  this 
procedure  is  termed  a  sequential  or  simplified  joint  inversion  process  in 
Chapter  2. 

For  TFM  and  IDEM  datasets,  a  proposed,  prototype  algorithm  for 
integrated  interpretation  is  illustrated  in  Figure  23.  The  algorithm  assumes 
that  the  TFM  data  is  converted  into  the  magnitude  of  the  analytic  signal 
(e.g.,  Hollyer,  Racic,  and  Butler  1997),  as  in  the  example  in  Figure  18.  The 
magnitude  of  the  analytic  signal,  equal  to  the  square  root  of  the  sum  of  the 
squares  of  the  three  spatial  gradients  of  the  TFM  (Roest,  Verhoef,  and 
Pilkington  1992),  is  determined  from  the  TFM  data  by  first  applying  a 
convolutional  filter  to  the  gridded  TFM  data  to  obtain  the  horizontal  gradients 
and  then  a  fast  Fourier  transform  or  a  discrete  Hilbert  transform  (Butler  1995) 
to  obtain  the  vertical  gradient  (Hollyer,  Racic  and  Butler  1997;  Geosoft  1996). 
The  advantage  of  the  TFM  analytic  signal  is  that  dipolar  anomalies  are 
converted  to  monopolar  anomalies  centered  over  the  buried  object,  making 
automatic  peak-picking  feasible  and  more  efficient.  Coincident  peaks  selected 
by  the  algorithm  will  be  related  to  ferrous  objects,  for  which  the  depths  are 
computed  by  analyses  of  the  TDEM  anomaly  characteristics  (Hollyer,  Racic, 
and  Butler  1997). 

Interpreted  depth  (to  center  of  the  ferrous  object,  Zo)  from  the  TDEM 
analysis  is  used  as  a  fixed  parameter  (constraint)  in  the  inversion  of  the 
magnetic  analytic  signal  anomaly.  The  value  of  using  the  TDEM  depth  to 
constrain  the  magnetic  inversion  is  that  the  fewer  the  number  of  “free* 
parameters  in  an  inversion,  the  faster  the  inversion  will  converge  and  the  more 
“unique*  the  inversion  result  (model)  will  be.  Alternatively,  the  prototype 
algorithm  can  be  modified  to  obtain  the  horizontal  location  (Xo,  Yo)  as  well  as 
the  depth  (Zo)  from  the  TDEM  analysis,  to  further  constrain  the  magnetic 
analytic  signal  inversion.  Finally,  inverting  the  magnetic  analytic  signal  gives 
the  magnetic  moment,  which  can  be  correlated  to  the  ferrous  mass  of  the  object 
causing  the  anomaly.  The  ferrous  mass  is  then  correlated  to  ferrous  masses  of 
known  ordnance  items  to  produce  a  likely  UXO  type  classification  (e.g., 
Pennella  1982).  This  prototype  algorithm  (Figure  23)  can  be  implemented  with 
minimal  effort  in  existing  data  analysis  and  visualization  software  packages 
(e.g.,  Geosoft  1996). 

A  modified  prototype  algorithm,  illustrated  in  Figure  24,  analyzes  the  TFM 
data  directly,  using  the  TDEM-determined  location  and  depth  as  constraints, 
and  determines  a  best-fitting  prolate  spheroid  model.  A  manual  implementation 
example  of  this  modified  prototype  algorithm  is  demonstrated  in  the  next 
chapter.  Variations  in  this  proposed  approach  to  joint  interpretation  are 
obvious.  For  example,  constraints  resulting  from  interpretation  of  TFM  data 
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can  be  applied  to  inversion  of  the  TDEM  data  to  a  final  model.  Also,  the 
procedure  can  be  iterated  to  minimize  an  overall  goodness  of  fit  criteria  to  the 
measured  data,  closely  approaching  a  true  joint  inversion. 
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Figure  23 .  Prototype  algorithm  for  integrated  multisensor  data  processing 
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Figure  24.  Modified  prototype  algorithm  for  integrated  multisensor  data  processing 


6  Development  of 
Phenomenological 
Modeling  Capability 


The  value  and  role  of  geophysical  signature  modeling  capabilities  are  dis¬ 
cussed  in  Chapter  2.  A  major  effort  of  this  study  was  development  of  analytical 
modeling  capability  for  geometries  that  more  closely  approximates  ordnance 
than  the  sphere.  Primary  weaknesses  of  the  spherical  model  of  ordnance  is  that 
the  spherical  synunetry  gives  no  orientation  dependence  of  the  ordnance  relative 
to  the  inducing  field.  Also,  for  magnetic  signature  modeling,  the  ordnance  is 
assumed  to  be  a  solid  ferrous  sphere,  with  the  size  determined  from  the  ferrous 
mass  and  the  density  of  iron  or  steel,  and  the  sphere  cannot  reflect  demag¬ 
netization  effects  associated  with  the  elongated  ordnance  geometry.  While  the 
geometry  of  ordnance  could  be  modeled  nearly  exactly  with  finite  element  or 
other  discretized  numerical  modeling  approaches,  the  computational  time  and 
the  parameters  of  the  model  may  not  be  conducive  to  inverse  modeling.  Dis¬ 
cretized  numerical  modeling  is  valuable  for  studying  phenomenological  effects, 
but  the  thrust  of  the  current  modeling  efforts  are  to  develop  analytical  modeling 
tools  which  can  be  adapted  ultimately  for  automated  inverse  modeling.  Simple 
geometries  which  more  closely  approximate  ordnance  shapes  than  the  sphere  are 
the  prolate  spheroid  and  circular  cylinders.  Both  full-field  analytical  and 
multipole  expansion  solutions  for  &e  induced  magnetic  field  of  a  prolate  spher¬ 
oid  (or  oblate  spheroid)  exist;  this  geometry  is  exploited  in  the  current  TFM 
modeling  capability  development.  TDEM  modeling  for  geometries  other  than 
the  sphere  is  difficult  and  is  not  developed  to  the  same  extent  as  TFM  modeling 
capability. 


Magnetic  Modeling  Program 


Background 

Altshuler  (1996)  reviews  the  deficiencies  of  a  spherical  model  for  ordnance 
and  evaluates  a  prolate  spheroid  model.  A  key  conclusion  of  Altshuler  from 
studies  of  solid  versus  spherical  shell  and  spheroidal  shell  models  is  that  it  is 
the  outer  ferrous  volume  of  the  ordnance  that  is  critical  to  the  magnetic 
signature,  not  the  ferrous  mass.  The  prolate  spheroid  model  is  a  realistic 
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representation  of  the  general  shape  of  ordnance  and  also  has  the  elongated 
geometry  of  ordnance  that  can  replicate  demagnetization  and  orientation  effects. 
A  prolate  spheroid  with  the  length  and  diameter  of  an  ordnance  item  is  a  good 
approximation  to  the  outer  ferrous  volume  of  the  ordnance.  The  induced  mag¬ 
netic  field  external  to  a  spheroid  can  be  determined  by  a  full  field  solution 
(Stratton  1941;  Altshuler  1996)  or  a  multipole  expansion  (McFee  and  Das 
1990).  In  the  multipole  expansion,  there  is  no  monopole  term,  and  the  quadra- 
pole  term  is  zero  due  to  symmetry.  Thus,  an  octapole  term  is  the  next  higher 
term  after  the  dipole  term;  and,  since  the  octapole  term  falls  off  as  1/r*,  there  is 
no  practical  reason  to  include  anything  higher  than  the  octapole  term.  Altshuler 
(1996)  compares  the  dipole  field  (for  a  spheroid  model)  with  the  full  field 
solution  and  concludes  that  for  measurement  distances  greater  than  about  two 
semimajor  axes  (the  length)  from  the  center  of  volume  of  the  spheroid,  the 
dipole  model  field  prediction  is  within  10  percent  of  the  full  field  model  predic¬ 
tion.  Thus  for  small  distances,  the  octapole  field  contribution  becomes  more 
significant.  Advantages  of  a  multipole  solution  compared  to  a  full  field  solution 
are  primarily  (a)  reduced  computational  time  and  (b)  ability  to  separate  the 
prolate  spheroid  dipole  term  for  comparison  with  the  dipole  solution  for  a 
sphere;  with  the  primary  disadvantage  being  a  possible  minor  lack  of  accuracy 
for  very  close  distances  of  the  model  to  the  signature  calculation  plane. 

The  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
implementation  of  a  multipole  expansion  magnetic  field  solution  for  a  spheroid 
model  follows  the  outline  of  McFee  and  Das  (1990).  Geometry  for  the 
magnetic  modeling  problem  is  illustrated  in  Figure  8.  There  are  two  obvious 
coordinate  systems  for  the  spheroid  model,  a  body-centered  system  (x„X2,X3), 
with  the  Xj  axis  along  the  long  symmetry  axis  of  Ae  spheroid,  and  a  space- 
centered  system(Xi  ,X2  ,X3 ),  in  which  the  spheroid  center  is  at  (Xo,Yo,2Io)  as 
shown  in  Figure  8.  Assuming  uniform,  parallel  ma^etization  M  of  the 
spheroid,  the  multipole  expansion  of  the  secondary  induced  magnetic  induction 
at  a  point  r  in  free  space  is 

b  =  b®  -I-  b®  ,  (1) 

where  b®  is  the  dipole  field  component  and  b®  is  the  octapole  field  component. 
The  dipole  field  component  is  given  by 

b®  =  [mo/  (4711^)]  { -m®  -I-  (3/r^  (r  •  m® )  r },  (2) 

where  no  is  the  magnetic  permeability  of  free  space,  r  =  |  r  |  ,  and  m®  is  the 
dipole  moment.  The  expression  for  the  octapole  field  b®  is  more  complicated, 
involving  a  rank  three  octapole  moment  tensor.  The  dipole  moment  is  related 
to  the  magnetization  by 

m®  =  MV  ,  (3) 

where  V  =  (4/3)Tiea®  is  the  volume  of  a  spheroid,  e  =  L/  (2a)  is  a  shape  factor 
for  the  spheroid,  L  is  the  length  of  the  spheroid,  and  2a  is  the  spheroid  dia¬ 
meter  orthogonal  to  the  symmetry  axis.  This  formulation  is  also  valid  for  a 
sphere,  where  e  =  1  (i.e.,  L  =  2a  and  a  is  the  radius  of  the  sphere)  and 
V  =  (4/3)tta^  is  the  volume  of  a  sphere.  Also,  for  the  sphere,  in  die  absence 
of  permanent  magnetization,  M  will  always  be  in  the  direaion  of  the  inducing 
earth’s  field  and  the  octapole  field  component  b®  is  zero. 
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For  the  spheroid,  the  induced  field  will  not  always  be  in  the  direction  of  the 
earth’s  inducing  field  even  in  the  absence  of  permanent  magnetization.  Resolv¬ 
ing  the  earth’s  field  bo  into  components  along  the  body-centered  coordinate 
system  (boi,bo2,bo3),  the  magnetization  in  the  spheroid  can  be  expressed  (Stratton 
1941,  McFee  and  Das  1990) 

M  =  (1/  Ho)  F  •  bo ,  (4) 

where  the  diagonal  components  F; ,  I  =  1,2,3,  of  the  demagnetization  tensor  F 
(the  tensor  is  diagonalized  in  the  body-centered  coordinate  system)  are  called 
the  demagnetization  factors.  Furthermore,  for  the  spheroid,  Fj  =  Fj .  The 
demagnetization  factors  are  functions  of  the  shape  factor  e  and  the  relative  mag¬ 
netic  permeabilities  of  the  spheroid  material  (Hti~  fh)  surrounding 

medium  ( Mrz  =  M2  /  Mo  )•  For  the  case  of  an  ordnance  item  buried  in  soil, 
generally  Hti  >  >  Mr2  and  M12  ®  1-  Finally,  for  a  sphere,  e  =  1  and  Fj  =  Fj 
=  Fj .  For  &e  spheroid,  two  cases  can  be  distinguished,  the  prolate  spheroid 
with  e  >  1  (the  case  used  to  model  UXO),  and  the  oblate  spheroid  with  e  <  1 
(which  might  be  used  to  model  a  ferrous  metallic  mine).  The  modeling  program 
checks  for  the  value  of  e  and  uses  the  appropriate  expressions  for  Fj .  Also,  the 
program  allows  for  calculation  of  b  (Equation  1)  or  just  the  dipole  field 
component  b®. 

Finally,  the  calculated  components  in  the  body-centered  coordinate  system 
are  rotated  into  the  space  or  earth-referenced  coordinate  system 

B  =  A  •  b  ,  (5) 

where  A  is  the  Euler  rotation  tensor  (e.g.,  Arfken  1985).  The  space  coordinate 
system  employed  for  the  solution  (and  the  computer  program)  is  a  right-handed 
system  with  X positive  to  the  north,  Y positive  to  the  -west,  and  Z positive  up¬ 
ward.  For  a  buried  ordnance  item  model,  the  depth  to  the  center  of  the  spher¬ 
oid  is  negative.  Inclination  of  the  spheroid  is  relative  to  the  positive  Z-  axis, 
and  thus  an  inclination  of  90  deg  is  horizontal.  Azimuth  of  the  spheroid  is 
defined  as  the  angle  between  a  vertical  plane  through  the  spheroid  and  the 
positive  X-axis  (north),  with  counterclockwise  angles  positive.  This  modeling 
convention  contrasts  to  the  customary  definitions  of  inclination  and  azimuth 
used  for  the  validation  and  phenomenology  studies:  inclination  (I)  positive  down 
from  horizontal;  azimuth  (A)  positive  counterclockwise  from  north. 

The  program  allows  calculation  and  plotting  of  the  magnitude  of  the  secon¬ 
dary  field  B  (Equation  5),  the  magnitude  of  the  total  field  B,  =  B  -I-  Bq,  the 
components  of  the  total  field  B,;  =  B;  -h  Ba  (I  =  1,2,3)  that  would  be  measured 
by  vector  magnetometers,  and  the  magnitude  of  the  total  magnetic  field 
anomaly,  i.e.,  the  value  measured  by  total  field  magnetometers  IBJ  minus  the 
ambient  Earth’s  field. 


ABtmf  =  I  B  -I-  Bo  I  -Bo  =  (B^  +  Bo'  -f  2  Bo*  B  -Bo,  (6) 

where  the  A  indicates  an  anomaly  value.  The  octapole  component  contribution 
to  the  preceding  values  calculated  by  the  magnetic  modeling  program  can  be 
“switched  off”  for  plotting  just  the  dipole  contribution. 
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Program  description  and  input/output 

The  magnetic  modeling  program  is  implemented  in  the  MATLAB*  (regis¬ 
tered  Trademark  of  The  Math  Works,  Inc.,  Natick,  MA)  programming  envi¬ 
ronment  to  take  advantage  of  the  built-in  graphics  capabilities  during  the 
development  and  validation  phases  of  the  capability.  The  program  is  initiated 
by  executing  the  script  magmod.  A  graphical  user’s  interface  (GUI)  simplifies 
and  expedites  input,  model  execution,  and  specification  of  graphical  output  for¬ 
mat.  Model  input  parameters  and  output  format  are  specified  on  three  GUI 
screens.  The  complete  program  listing  is  at  Appendix  B. 

Main  screen.  The  first  screen  is  the  main  or  master  screen  that  has  several 
“push  buttons”  for  accomplishing  the  following: 


Push  Button 

Description 

Edit  input  Parameters 

Enters  the  second  screen,  for  data  input  specifications 

Run  Mode! 

Executes  the  model  calculations,  with  the  current  input  parameters 
and  plot  specifications 

Edit  Piot  Instructions 

Enters  the  third  screen  to  specify  output  plot  format,  including  page 
orientation,  number  of  plots,  types  of  plots,  etc. 

Plot  Results 

Reproduces  the  specified  graphical  output  on  the  computer  screen 

Hard  Copy 

Produces  a  hard  copy  of  the  output  on  the  specified  printer/pl otter 

In  addition,  there  are  three  secondary  push  buttons  for  specifying  a  “batch  run- 
plot,”  to  make  a  movie  from  the  output  graphics,  and  to  play  the  movie.  Two 
fill-in  boxes  allow  specification  of  a  comparison  data  file  and  a  code  for  the 
desired  printer  or  plotter.  The  comparison  data  file  could  be  TFM  profile  mea¬ 
surements  acquired  over  known  or  unknown  anomaly  targets  with  the  same 
distance  scale  as  the  calculated  profile  values.  The  comparison  data  profile  can 
be  superimposed  on  the  calculated  profile;  this  is  a  useful  feature  for  program 
validation  or  for  iterative  forward  modeling. 

Input  parameter  screen.  The  input  parameter  screen  consists  of  a  number 
of  fill-in  boxes  that  can  be  freely  edited.  To  the  right  side  of  each  fill-in  box  is 
a  description  of  the  parameter  to  be  entered.  Initially  there  is  a  check  box  to 
indicate  metric  units  (or  English  units  if  unchecked).  The  following  tabulation 
gives  the  parameter  list  in  the  order  of  appearance  on  the  screen: 

Input  Parameter  Description 

_  Minimum  X-  (north)  value  of  grid  in  m  or  ft 

_  Maximum  X-  (north)  value  of  grid  in  m  or  ft 

_  Minimum  Y-  (west)  value  of  grid  in  m  or  ft 

_  Maximum  Y-  (west)  value  of  grid  in  m  or  ft 

_  Height  of  instrument  (sensor)  above  ground 

in  m  or  ft 

_  Earth’s  magnetic  field  in  nanoTeslas  (nl) 

_  Inclination  of  earth’s  field  in  degrees 

*.mat  Output  file  name 

_  Length  of  body  in  m  or  in. 
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_  Diameter  of  body  in  m  or  in. 

_  Relative  magnetic  permeability  of  body 

_  Angle  of  body  symmetry  axis  relative  to  Z-axis 

in  degrees 

_  Angle  of  body  vertical  plane  relative  to  X  in 

degrees,  positive  counterclockwise  (CCW) 

_  X-coordinate  QQ  of  body  center  in  m  or  ft 

_  Y-coordinate  (Yq)  of  body  center  in  m  or  ft 

_  Z-component  (Zo)  of  body  center  in  m  or  ft 

(negative  for  buried  objects) 

For  a  given  set  of  anomaly  calculations  for  a  list  of  models,  many  of  the  input 
parameters  will  remain  constant  between  program  runs. 

Plot  speciftcation  screen.  The  plot  screen  allows  specification  of  the 
numerous  plot  options  available.  At  the  top  are  two  check  boxes,  to  speciiy 
dipole,  octapole,  or  both  components.  The  Plot  Format  click  box  brings  up  a 
number  of  portrait  and  landscape  plot  orientation  options.  The  File  Name  fill- 
in  box  is  used  to  specify  the  output  file  name  (same  as  specified  on  the  param¬ 
eter  input  screen);  this  box  is  filled  in  by  default  with  the  most  recently 
completed  model  calculation  result  but  may  be  edited  to  plot  other  saved  output 
files.  As  many  as  six  plots  or  graphs  can  be  specified  for  a  single  output  page; 
these  are  indicated  beneath  the  heading  “Data  Type.”  Clicking  on  any  of  the 
six  boxes  brings  up  a  list  of  data  types  which  can  be  plotted,  and  highlighting  a 
parameter  selects  it  for  plotting.  The  data  types  which  can  be  plotted  are: 

Data  Type  Plot  Type  Cross  Yalaes 

~  X-component  of  total  magnetic  field 

~  Y-component  of  total  magnetic  field 

~  Z-component  of  total  magnetic  field 

-  Magnitude  of  the  secondary  field 

-  Total  field  anomaly 

-  Vector  representation  of  horizontal  components 

-  Representation  of  spheroidal  body 

~  Surface  of  spheroid  body 

~  Listing  of  major  parameters 

For  the  data  types,  except  the  last  three  options  in  the  above  list,  selection 
brings  up  a  click  box  under  the  heading  Plot  Type. 

Clicking  a  Plot  Type  box  opposite  a  selected  Data  Type  brings  up  a  list  of 
options  for  how  to  plot  the  data.  Various  types  of  plots  are  possible;  3-D 
surface  representations;  contour  plots;  profile  plots  (projections).  Three- 
dimensional  surface  plots  are  representations  of  the  value  of  a  calculated 
magnetic  value  (along  a  vertical  axis)  and  its  variation  over  the  specified  ranges 
of  X  and  Y.  The  3-D  plots  can  be  color-filled,  a  shade  of  gray,  shaded  relief, 
or  have  superimposed  contours.  Likewise,  the  contour  plots,  which  are  two- 
dimensional  (2-D)  can  be  a  shade  of  gray,  have  color  contour  lines,  or  color- 
filled  contours.  An  effective  display  option  is  to  have  a  3-D  representation  plot 
vertically  offset  from  its  2-D  contour  plot  counterpart.  The  projection  or 
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profile  plot  option  is  a  familiar  way  of  visualizing  the  detailed  variation  and 
shape  of  the  anomaly  signature  along  a  surface  line  parallel  to  the  X-axis  (north- 
south)  or  parallel  to  the  Y-axis  (east-west).  From  one  to  four  profile  or  pro¬ 
jection  lines  can  be  plotted  on  the  same  graph.  When  the  projection  plot  type 
option  is  seleaed,  a  box  opens  under  Cross  Values  to  specify  the  desired  loca¬ 
tions  of  the  profile  lines.  For  example,  if  the  Plot  Type  selection  “3  Projections 
Parallel  to  X-Axis”  is  selected,  three  Y-values  must  be  provided  in  the  Cross 
Values  box,  indicating  Y-values  where  the  profiles  cross  the  Y-axis.  A  common 
procedure,  for  the  above  example,  is  to  have  the  center  of  the  spheroid  beneath 
the  point  (0,0)  and  to  have  one  profile  cross  at  Y  =  0  and  the  remaining  two 
profiles  cross  at  Y  =  -1  m  and  Y  =  1  m.  This  procedure  gives  a  central  profile 
that  generally  indicates  the  maximum  signature  magnitudes  and  the  spatial 
wavelength  of  the  signature,  and  then  indicates  changes  in  signature 
characteristics  laterally;  the  profiles  allow  rational  consideration  of  along-track 
and  cross-track  measurement  spacings  required  to  capture  the  relevant  charac¬ 
teristics  of  the  signature  for  the  specific  case. 

Example.  A  specific  example  illustrates  most  of  the  input  and  output  fea¬ 
tures  discussed  above.  Results  of  a  model  validation  comparison  to  measure¬ 
ment  data  from  the  NRL  ordnance  signature  library  is  shown  in  Figure  25.  The 
comparison  is  for  a  105-mm  artillery  projectile,  oriented  horizontally  and 
pointing  north.  The  figure  consists  of  four  plots  corresponding  to  various  Data 
Type  and  Plot  Type  specifications:  (a)  the  upper  left  display  is  a  shaded  relief 
3-D  TFM  anomaly  plot  vertically  above  a  standard  contour  plot;  (b)  the  upper 
right  display  is  a  representation  of  the  input  spheroidal  model;  (c)  the  lower  left 
display  is  a  north-south  profile  TMF  anomaly  plot  (parallel  to  the  north-south 
axis  and  crossing  the  east-west  axis  at  0)  showing  model  calculation  (solid  line) 
and  superimposed  measurement  results  (+);  (d)  the  lower  right  display  is  a 
listing  of  the  program  input  parameters.  Also  included  in  the  lower  left  plot  is  a 
legend  displaying  the  profile  axis-crossing  value,  the  maximum  and  minimum 
TFM  values  along  the  profile,  and  the  profile  distance  between  maximum  and 
minimum  values. 


Magnetic  modeling  program  validation 

Validation  of  the  WES  magnetic  modeling  program  addresses  two  consider¬ 
ations,  reality  checks  and  correlation  or  comparison  with  well-controlled 
magnetic  signature  measurements. 

a.  Reality  checks. 

(1)  For  given  input  parameters,  the  calculated  anomaly  signature  must 
become  smaller  in  magnimde  and  broader  in  spatial  extent  as  the 
depth  of  the  spheroid  is  increased. 

(2)  For  a  spheroid  model,  identical  results  should  be  obtained  for 
indistinguishable  cases  with  different  input  parameter  values. 

(3)  Calculations  for  pairs  of  indistinguishable  inclinations,  e.g.,  0  and 
180  deg  or  90  and  270  deg,  must  be  identical. 

(4)  Calculations  for  pairs  of  indistinguishable  azimuths,  e.g.,  45  and 
225  deg,  must  be  indistinguishable. 
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Total  Field  Anomaly  Spheroid 
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Figure  25.  Example  output  of  magnetic  modeling  program.  The  total  field  anomaly  profile  plot  (lower  left)  includes  a  validation  model  calculation 

(solid  line)  of  a  specific  case  fi*om  the  NRL  ordnance  signature  library  (-f):  105-mm  artillery  projectile  oriented  horizontally  and  pointing  north 


(5)  For  a  prolate  spheroid  model,  the  calculations  must  show  changes  as 
the  object  is  rotated  (i.e.,  the  inclination  and/or  azimuth  changes) 
relative  to  the  Earth’s  field  orientation,  refiecting  demagnetization 
effects  due  to  the  elongated  geometry. 

(6)  For  the  simple  case  of  a  sphere  (i.e.,  L  =  2a  or  e  =  1),  the 
calculated  magnetic  anomaly  should  have  a  simple  dipolar  shape  and 
should  not  vary  as  the  input  angles  are  varied. 

b.  Correlation  with  well-controlled  magnetic  signature  measurements. 

An  example  of  this  type  of  validation  is  given  in  Figure  25,  where 
the  calculated  signature  for  a  spheroidal  model  of  a  105-mm  projectile 
agrees  very  closely  with  measured  signature  data  over  an  inert  105-mm 
projectile.  This  type  of  detailed  profile  comparison  was  conducted  for 
numerous  cases  from  the  NRL  ordnance  signature  database.  Good 
correlation  of  calculated  signatures  with  all  cases  from  the  signature  database 
is  achieved  with  the  specified  length  and  diameter  of  the  ordnance  as  input 
parameters  for  the  spheroid,  except  for  the  60-  and  81-mm  mortars,  where 
the  length  is  modified  to  reflect  the  ferrous  length  and  not  the  total  ordnance 
item  length.  Detailed  comparisons  of  magnetic  modeling  program  signature 
calculations  with  measured  signatures  from  the  WES  ordnance  signatures 
database  are  also  good  and  will  be  presented  in  a  subsequent  report. 

Another  validation  example  is  shown  in  Figure  26,  where  color-filled 
contour  maps  of  calculated  signatures  (bottom)  are  compared  directly  with 
measured  signatures  (top)  for  a  105-mm  projectile  at  five  azimuths  in  the 
horizontal  plane.  The  agreement  between  the  top  and  bottom  in  both 
specifics  and  generalities  is  quite  good.  Note  that  the  calculated  signatures 
for  azimuths  of  0  and  180  deg  are  identical  as  expected.  The  measured 
signatures  for  the  0-  and  180-deg  azimuth  cases  do  show  some  differences, 
however  it  is  not  possible  to  determine  from  this  one  example  if  the 
differences  are  due  to  the  fore/aft  asymmetry  of  the  actual  ordnance  item  or 
is  due  to  noise  in  the  180-deg  azimuth  measured  case. 


Phenomenology  studies  with  magnetic  modeling  program 

After  validation,  the  magnetic  modeling  program  can  be  used  with  some 
confidence  for  phenomenology  studies  and  for  multisensor  integration  investi¬ 
gations.  A  valuable  use  of  a  magnetic  modeling  program  is  to  study  or  predict 
the  detectability  of  a  given  ordnance  item  as  a  function  of  depth  and  orientation 
and  to  compare  the  detectability  of  different  ordnance  items.  For  example, 
graphs  in  Figures  27  and  28  compare  the  maximum  positive  TFM  anomaly  for  a 
155-mm  projectile  and  a  500-lb  bomb,  respectively,  as  a  function  of  depth  to 
object  center  for  three  inclination  and  azimuth  combinations  (Note  the 
definitions  for  inclination  (I)  and  azimuth  (A)  in  Figures  27  and  28  used  for  all 
validation  and  phenomenology  studies  in  this  report).  With  an  estimate  of  site- 
specific  magnetic  noise,  which  includes  both  time-varying  and  background 
“clutter”  components  (including  spatial  variability),  maximum  depths  for 
detection  of  given  ordnance  items  can  be  estimated.  For  example,  if  the 
average  magnetic  noise  level  at  a  site  is  5  nT  and  assuming  that  measurement 
spacing  is  adequate  to  delineate  the  anomaly  signature,  then  a  10-nT  anomaly 
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Figure  26.  Validation  example  for  measured  (top)  and  calculated  (bottom)  TFM  signatures  for  105-mm  artillery  projectile,  oriented  horizontally 
and  for  azimuths  of  0, 45,  90,  1 35, 1 80  deg.  Depth  below  surface  is  0.5  m  and  instrument  height  is  0.25  m 
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Figure  27.  Maximum  positive  TFM  anomaly  calculations  for  1 55-mm  artillery  projectile  for  three  inclination  and  azimuth  orientation  combinations 
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Figure  28.  Maximum  positive  TFM  anomaly  calculations  for  500-lb  bomb  for  three  inclination  and  azimuth  orientation  combinations 


(two  times  the  average  noise  level)  is  reliably  detectable.  Using  the  10-nT 
threshold  for  detection  and  Figure  27,  a  155-mm  projectile  is  detectable  to  a 
depth  of  approximately  2  m  for  the  worst-case  orientation  (inclination  =  0, 
relative  to  horizontal,  and  azimuth  =  90  deg)  and  to  a  depth  of  approximately 
4  m  for  the  best-case  orientation  (inclination  =:  90  deg).  Similarly,  using  the 
10-nT  threshold  and  Figure  28,  a  500-lb  bomb  is  detectable  to  a  maximum 
depth  of  approximately  4  m  for  the  worst-case  orientation  and  to  a  maximum 
depth  of  approximately  9  m  for  the  best-case  orientation. 

Another  example  of  phenomenology  studies  is  illustrated  in  Figure  29.  This 
example  dramatically  illustrates  the  effects  of  orientation/demagnetization 
effects  caused  by  the  orientation  of  the  elongated  ordnance  (155-mm  projectile) 
relative  to  the  earth’s  magnetic  field.  For  both  cases  in  Figure  29,  the  inclina¬ 
tion  of  the  earth’s  field  is  65  deg,  which  is  typical  of  the  inclination  at  Fort 
Carson,  CO,  and  Fort  A.  P.  Hill,  VA.  In  the  first  case  (left),  the  inclination  of 
the  spheroid  model  (45  deg  relative  to  horizontal)  is  in  the  approximate 
direction  of  the  earth’s  field,  and  the  TMF  anomaly  signature  is  dipolar  in 
appearance  with  a  peak  positive  magnitude  of  nearly  240  nT.  For  the  second 
case  (right),  the  inclination  of  the  spheroid  model  (-45  deg  relative  to  horizon¬ 
tal)  is  approximately  perpendicular  to  the  Earth’s  field,  and  the  TFM  anomaly 
signature  is  monopolar  in  appearance  with  a  peak  positive  magnitude  of  only 
110  nT  (less  than  one-half  the  magnitude  of  the  first  case).  This  example  makes 
it  dramatically  obvious  why  a  simple  magnetic  dipole  (equivalent  sphere) 
inversion  can  result  in  significant  errors  in  ordnance  depth  and/or  size. 


Electromagnetic  Induction 
Modeiing  Considerations 


Background 

Due  to  the  complexity  of  general  solutions  of  EMI  responses  for  realistic 
UXO  geometries,  this  section  only  briefly  reviews  work  of  others  and  pre¬ 
liminary  investigations  conducted  during  the  present  effort.  Rigorous  analytical 
solutions  to  date  are  for  simple  geometries,  e.g.,  layered  earth,  spheres,  plates, 
and  circular  cylinders  (McNeill  and  Bosnar  1996,  Barrow  et  al.  1996,  Ward 
and  Hohmann  1987).  Das  et  al.  (1990)  extend  consideration  to  a  prolate 
spheroid  but  point  out  that  rigorous  analytical  solutions  do  not  exist  for  the 
general  case  of  a  conducting,  permeable  spheroid.  Analytical  solutions  often 
require  numerical  evaluation  or  approximation  but  do  not  require  problem 
discretization,  as  in  finite  difference,  finite  element,  boundary  integral,  or 
hybrid  computational  techniques.  Numerical  modeling,  involving  space,  time 
and/or  frequency  discretization,  complex/realistic  geometries,  and  the  full 
range  of  physical  properties  and  processes,  is  notoriously  computationally 
intensive  (Hohmann  1987,  Laveley  1996);  currently  this  type  numerical 
modeling  is  reserved  for  phenomenological  studies  but  is  not  practical  for  the 
repetitive,  iterative  calculations  generally  required  for  geophysical  inversion. 
Thus  numerical  modeling  and  even  analytical  solutions  resort  to  approxima¬ 
tions  of  various  types:  asymptotic  approximations;  quasi-static  solutions;  low- 
frequency  or  high-frequency  limit  solutions  for  FDEM;  early-time  or  late-time 
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Figure  29.  Example  of  phenemonology  study.  Magnetic  model  signature  calculation  for  155-mm  artillery  projectile  with  inclination  (a)  approximately  in  the 
direction  of  the  earth’s  field  inclination  and  (b)  approximately  perpendicular  to  the  earth’s  field 


solutions  for  TDEM.  Many  of  these  approximations  in  time  or  frequency 
domains  are  equivalent  for  different  combinations  of  physical  properties  (Wait 
1982,  Butler  and  Fitterman  1986,  Butler  1986,  Ward  and  Hohmann  1987). 

As  discussed  in  Chapter  2,  EMI  methods,  particularly  TDEM  methods,  are 
used  extensively  for  UXO  detection  and  mapping  surveys.  An  extensively  used 
TDEM  system,  the  Geonics  EM61  standard  system  consists  of  a  1-m  x  1-m 
transmitter  (Tx)  coil  and  two  1-m  x  1-m  receiver  (Rx)  coils,  with  one  Rx  con¬ 
centric  to  the  Tx  and  the  second  Rx  offset  0.4  m  vertically  (Figure  30a).  A  spe¬ 
cialized,  nonstandard  0.5-m  x  0.5-m  Tx  and  Rx  coil  EM61  system  was  used  in 
the  DARPA  program.  The  EM61,  like  many  TDEM  systems,  operates  by 
rapidly  turning  off  a  current  in  the  Tx  and  detecting  transient  (decaying)  magne¬ 
tic  fields  at  the  Rx’s  from  induced  current  transients  in  the  subsurface.  The 
physics  of  the  transient  induction  process  and  descriptions  of  general  purpose 
TDEM  systems  are  reviewed  by  Butler  and  Fitterman  (1986),  and  the  induction 
process  in  compact,  conductive,  permeable  metallic  objects  outlined  descrip¬ 
tively  by  McNeill  and  Bosnar  (1996).  A  schematic  of  the  transmit,  indurtion, 
and  measurement  process  is  illustrated  in  Figure  30b-d. 

General  purpose  TDEM  systems  sample  the  transient  signal  with  many  time- 
gates  or  channels,  typically  20  to  30.  The  standard  EM61  measures  the  vertical 
component  of  the  secondary  field  and  integrates  the  transient  over  one  time-gate 
extending  from  0.47  to  0.87  ms;  the  time  gate  is  selected  to  enhance  sensitivity 
to  metallic  objects  (for  the  NRL  MTADS  system  described  previously,  the  time- 
gate  is  moved  closer  to  the  transmit  pulse  for  even  higher  sensitivity;  Nelson, 
McDonald,  and  Robertson,  1997).  The  prototype  TDEM  system  (EM61-3D) 
used  at  the  DARPA  sites  records  20  time  channels  for  each  of  3  orthogonal  Rx 
coils  (the  vertical  and  2  horizontal  components). 

A  decay  curve  recorded  with  the  vertical  component  Rx  coil  of  the  prototype 
EM61-3D  system  is  shown  in  Figure  31,  with  the  time-gate  (window)  of  the 
standard  EM61  indicated.  The  decay  curve  in  Figure  31  was  recorded  over  a 
20-cm  X  20-cm  aluminum  plate  buried  at  5  cm  (approximately  45  cm  below  the 
Tx)  at  one  of  the  DARPA  sites.  The  two  Rx  measurements  for  the  standard 
EM61  give  two  integrated  magnitude  measurements  and  their  difference  (or 
vertical  gradient)  as  a  function  of  position  as  the  system  is  pulled  along  the 
ground  surface.  Typically,  the  EM61  data  are  displayed  in  map  form,  and 
signatures  of  shallow  metallic  objects  are  apparent,  allowing  accurate  location 
of  the  positions  of  the  objects.  Depth  estimates  are  determined  from  the  spatial 
wavelengths  of  the  signatures  or  based  on  gradient  information  from  the  top  and 
bottom  Rx’s.  Current  efforts  apply  empirical  adjustments  based  on  measured 
signatures  to  spherical  model  predictions  to  account  for  ordnance  shape  and 
orientation  effects  (Barrow  et  al.  1996  and  1997)  with  increasing  success. 

For  the  limited  information  acquired  with  the  standard  EM61  and  typical 
FDEM  systems  used  for  UXO  surveys  that  operate  at  only  one  or  two  fre¬ 
quencies,  detection,  location,  depth  estimates,  and  limited  discrimination  capa¬ 
bility  likely  represent  the  maximum  capability.  Additional  capability  for 
discrimination  is  possible  when  the  EMI  data  are  integrated  with  other  geophys¬ 
ical  data  (Chapter  2).  Another  possibility,  for  increasing  discrimination  and 
identification/classification  capability,  is  to  obtain  more  information  with  the 
EMI  systems.  Two  enhanced  capability  prototype  EMI  systems  (first  used  for 
UXO  applications  during  the  DARPA  program)  are  described  in  more  detail  in 
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a.  Photograph  of  standard  EM61. 


b.  Current  in  Tx  coil 


c.  Induced  electromotive  force  caused 
by  Tx  current  turn-off 


d.  Secondary  (transient)  magnetic 
field  caused  by  eddy  currents 


Figure  30.  (a)  Photograph  of  the  standard  EM6 1 ;  (E-d)  schematic  of  the  transient  induction  and  measurement 
process 
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Amplitude  (mv) 


Figure  3 1 .  Vertical  component  decay  curve  recorded  with  EM6 1  -3D  over  an  aluminum  plate  (20  x  20  cm) 
buried  at  5  cm  at  a  DARPA  site;  time  window  for  standard  EM6 1  is  indicated 
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the  next  chapter.  These  enhanced  capability  systems  include  a  multichannel, 
multicomponent  TDEM  system  and  a  multifirequency  FDEM  system.  EMI 
modeling  and  analysis  capability  during  the  current  effort  is  addressed  specifi¬ 
cally  to  the  multichannel,  multicomponent  TDEM  system;  since  it  likely  has  the 
most  general  applicability,  and  then  modeling  the  standard  EM61  signatures  of 
ordnance  is  a  special  case. 

Three  techniques  for  modeling  the  full  time  domain  EMI  response  of  a  com¬ 
pact,  conductive,  permeable  ordnance  object  are  (1)  simple  model-based 
approach  with  empirical  adjustments  for  time  response  (quasi-empirical 
approach),  (2)  complete  physics-based  analytical  time  domain  solution  for  well- 
defined,  realistic  geometry  (e.g.,  a  prolate  spheroid),  and  (3)  a  discretized, 
numerical  complete  solution  for  actual  ordnance  geometry.  Preliminary  inves¬ 
tigations  of  a  quasi-empirical  approach  are  presented  here.  The  capability 
allowed  by  approach  (2)  would  greatly  enhance  the  feasibility  of  true  Joint 
inversion  for  multisensor  integration  and  should  be  aggressively  pursued. 
Recently  initiated  efforts  address  approach  (3)  under  a  variety  of  programs; 
these  efforts  will  significantly  advance  understanding  of  TDEM  phenomenology 
and  enable  parameter  studies,  although  utility  of  these  approaches  as  tools  for 
data  inversion  is  uncertain. 


Quasi-empirical  TDEM  modeling  considerations 

McNeill  and  Bosnar  (1996)  discuss  the  response  of  compact,  simple  geom¬ 
etry  objects  to  static  and  transient  EM  fields.  Specifically,  they  reference 
spheres,  plates  (discs),  and  rods  (or  cylindrical  shells).  Also,  Ae  contributions 
of  finite  electrical  conductivity  and  finite  magnetic  permeability  to  the  overall 
EMI  response  are  considered.  To  model  these  simple  geometry  objects, 
appropriately  placed  and  oriented  magnetic  dipoles  reproduce  the  key  features 
of  the  quasi-static  spatial  signatures.  For  elongated  or  flattened  objects  such  as 
rods  and  discs,  the  EMI  response  is  distinctly  different  depending  on  the  orien¬ 
tation  of  the  object  relative  to  the  inducing  magnetic  field.  For  the  long  axis  of 
rods  perpendicular  to  the  inducing  field  or  the  plane  of  discs  perpendicular  to 
the  inducing  field,  the  response  of  permeable  (ferrous  material)  and 
nonpermeable  objects  are  similar;  when  the  inducing  field  is  abruptly  turned 
off,  a  dipolar  response  due  to  decaying  eddy  currents  (conductivity  response) 
dominates.  For  an  inducing  field  parallel  to  the  long  axis  of  rods  or  parallel  to 
the  plane  of  discs,  the  eddy  current  response  is  small,  but  a  dipolar-like 
polarization  response  is  induced  in  permeable  objects;  for  this  case  the  response 
is  virtually  zero  for  nonpermeable  objects.  For  a  general  orientation  of  the 
inducing  field,  the  response  of  plates  and  rods  can  be  approximated  by 
orthogonal  dipolar  responses  in  the  two  symmetry  directions  of  the  objects. 
McNeill  and  Bosnar  (1996)  present  measured  time  decays  of  the  induced  fields 
in  ferrous  and  non-ferrous  spheres,  plates,  and  cylindrical  shells  and  actual 
ordnance  items;  and  for  the  cylindrical  shells  and  ordnance  items,  the  responses 
for  the  inducing  field  parallel  to  and  perpendicular  to  the  long  axes  are  given. 
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McNeill^  proposes  an  initial  approach  for  time-domain  response  of  UXO  that 
utilizes  orthogonal  dipoles,  along  and  perpendicular  to  the  long  axis,  to  rep¬ 
resent  the  total  EMI  spatial  response  (see  also  Das  et  al.  1990).  The  procedure 
then  modifies  the  relative  contribution  of  each  dipole  according  to  the  measured 
decay  characteristics  of  specific  ordnance  items,  determined  with  the  ordnance 
item  parallel  to  and  perpendicular  to  the  primary  inducing  field  of  an  EM61  Tx. 
For  each  calculation  location  (X,Y),  the  inducing  field  (which  is  assumed 
uniform  at  the  object  position)  is  resolved  into  components  along  the  orthogonal 
dipole  orientations,  which  simulate  the  azimuth  and  inclination  of  the  ordnance. 
Indication  of  the  capability  of  this  approach  to  replicate  general  observations  of 
TDEM  signatures  and  spatial  signature  evolution  with  time  is  illustrated  in 
Figure  32  fi-om  McNeill*  (155-mm  projectile  simulated  at  depth  of  2  m, 
inclination  45  deg,  and  azimuth  45  deg).  The  upper  set  of  contour  plots  (2-m  x 
2-m  area)  is  the  root  square  magnitude  of  the  horizontal  field,  i.e.,  Bj,  =  (B/  -f- 
By^  )'^.  Both  the  horizontal  and  vertical  components  indicate  the  target  azimuth 
(45  deg).  The  effect  of  excitation  of  the  perpendicular  dipole  is  noted  at  early 
time,  which  decreases  with  time.  The  long  axis  excitation  dominates  at  late 
time. 

The  WES  implementation  of  this  quasi-empirical  approach  to  TDEM  signa¬ 
ture  modeling,  allows  further  investigation  of  the  signatures  as  funetions  of  the 
depth,  orientation  parameters,  and  time  evolution.  In  Figure  33  for  example, 
the  early  time  (0.3  ms)  dip  response  of  the  model  is  investigated,  for  a  depth  of 
2  m  and  azimuth  of  45  deg.  The  model  signatures  exhibit  an  intuitive  transition 
from  bi-orthogonally  symmetric  (symmetric  about  both  the  azimuth  direction 
and  perpendicular  to  the  azimuth)  for  a  dip  of  0  deg,  proceeding  to  axially 
symmetric  about  the  45  deg  azimuth  for  a  dip  of  45  deg,  and  finally  to 
completely  circularly  symmetric  for  a  dip  of  90  deg.  The  azimuth  response  of 
the  model  is  illustrated  in  Figure  34,  for  azimuths  of  45  and  90  deg  and  a  dip  of 
45  deg.  As  the  azimuth  is  changed,  the  signature  components  maintain  shape 
and  magnitude,  as  required.  Future  plans  for  this  model  include  efforts  to 
determine  rational  procedures  for  including  physical  size  of  the  object  in  the 
model  (currently  the  dipole  size  parameters  are  normalized  to  unity  and  size 
only  enters  in  the  measured  time  decay  over  actual  size  ordnance)  and 
validating  depth  and  orientation  signature  dependence  by  correlation  with 
measured  signatures. 


Phenomenological  Modeling 
and  Muitisensor  Data  Integration 


Geophysical  parameter  space  plots 

Selected  examples  of  multisensor  data  presented  previously  illustrate  data 
integration  procedures,  and  a  specific  type  of  analysis  in  geophysical 


‘  Personal  Communication,  July  1997,  J.  D.  McNeill,  Geonics,  Limited,  Mississauga,  Ontario, 
Canada. 
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Figure  32.  Theoretical  calculations  of  the  spatial  and  time  EMI  response  of  a  1 55-mm  ordnance  item.  Depth  =  2  m;  Azimuth  =  45  deg; 

Inclination  =  45  deg.  The  top  plots  are  the  root  mean  square  of  the  horizontal  components.  Each  plot  area  is  2  m  x  2  m  (McNeilF) 


I 


Figure  33.  Theoretical  early-time  (0.3  msec)  EMI  response  for  155-mm  projectile  model, 
depth  =  2  m,  azimuth  =  45  deg ;  model  dip  response 
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Time  0.3  msec 
Dip  45  degrees 


Horizontal  -  Components  -  Vertical 


Azimuth  ~  45  degrees 


Azimuth  =  90  degrees 

4-m  X  4-m  Plot  Areas 


Figure  34.  Theoretical  early-time  (0.3  msec)  EMI  response  for  155-mm  projectile  model, 
depth  =  2  m,  dip  (inclination)  =  45  deg;  model  azimuth  response 
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measurement/parameter  space  now  illustrates  the  importance  of  developing 
phenomenological  modeling  capability.  The  role  of  model-based  analyses  and 
interpretation  of  multisensor  geophysical  data  in  general  is  discussed  extensively 
in  Chapter  2.  An  example  of  TFM  and  TDEM  data  from  the  NRL  ordnance 
signature  library  in  Figure  35  serves  to  illustrate  the  key  points  (4  m  x  4  m 
area  for  top  plots).  The  data  are  for  a  105-mm  projectile  at  depth  0.49  m, 
inclination  0,  and  azimuth  90  deg.  The  TFM  data  in  the  left  plot  is  the  same  as 
one  of  the  cases  in  Figure  26,  used  to  validate  the  magnetic  modeling  program. 
The  TDEM  data  in  Figure  35  is  similar  and  consistent  with  the  TDEM  early- 
time  modeling  examples  in  Figures  32  through  34  (except  for  the  obvious 
difference  in  ordnance  size  and  inclination).  Thus  with  some  advances  in 
understanding  and  modeling  capability  for  the  TDEM  case,  both  of  the  sensor 
datasets  in  Figure  35  can  be  modeled,  allowing  a  model-based  interpretation. 
The  middle  plots  are  principal  N-S  profile  plots  across  the  center  of  the  signa¬ 
ture  in  the  top  plots.  The  bottom  plot  is  a  measured  parameter  space  plot, 
which  simultaneously  portrays  (integrates)  the  amplitude  and  phase  relationships 
of  the  two  signature  profiles.  Parameter  space  plots  for  four  cases  are  given  in 
Figure  36  for  two  orientations  of  the  same  ordnance  item  and  for  two  ordnance 
items  with  the  same  orientation;  depth  is  approximately  the  same  for  all  four 
cases. 

Parameter  space  plots  have  diagnostic  characteristics,  allowing  interpretation 
of  the  object(s)  causing  the  multisensor  anomalies  with  fewer  known  parameters 
and/or  assumptions  (e.g.,  Butler  1995).  Parameter  space  analyses  techniques 
are  not  limited  to  2-D  spaces,  with  a  3-D  space  from  three  datasets,  e.g.,  TFM, 
TFG,  and  TDEM,  a  common  possibility.  Also,  parameter  space  plots  can  be  in 
terms  of  extracted  key  parameters,  such  as  the  example  in  Figure  9.  In  terms  of 
multisensor  integration  by  nonmodel  based  approaches  such  as  neural  networks 
or  “fuzzy  logic,”  the  parameter  space  plot  collapses  multisensor  signature 
information  to  a  form  better  adapted  to  a  holistic  type  analysis.  For  model- 
based  approaches  to  multisensor  (multimethod)  integration,  assumption  of  a 
common  geometry  for  each  dataset  allows  direct  solution  for  the  common  model 
parameters,  e.g.,  length,  diameter,  depth,  inclination,  azimuth.  The  model- 
based  solution  proceeds  from  either  (a)  a  geometrical  analysis  of  the  parameter 
space  figure  or  (b)  iterative  forward  modeling  (inverse  solution  in  parameter 
space).  Forward  modeling  capability  for  each  geophysical  method,  once 
validated  as  discussed  previously,  allows  systematic  study  of  parameter  space 
geometries  in  a  more  rigorous  and  comprehensive  manner  than  relying  on 
measured  data  alone. 


Manual  implementation  of  proposed 
multisensor  integration  algorithm 

The  NRL  Twentynine  Palms  MTADS  data  afford  opportunities  for 
multisensor  integration  investigations.  Figure  37  shows  an  approximately  36-m 
X  36-m  area  display  of  TFM,  TFG,  and  TDEM  datasets  (McDonald  et  al. 

1997).  The  yellow  circles  in  Figure  37  indicate  the  location  of  a  group  of 
closely  spaced  inert  60-mm  mortar  rounds.  To  illustrate  the  role  of  phenomeno¬ 
logical  modeling  in  UXO  detection  and  discrimination,  a  manually  implemented 
version  of  the  proposed  algorithm  in  Figure  23  is  applied  to  the  geophysical 
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Figure  35.  Multisensor  datasets  from  the  NRL  ordnance  signature  library,  TFM  to  left  and 

IDEM  to  right;  105-mm  projectile  at  0.5  m  depth,  0  deg  inclination,  0  deg  azimuth. 
Lower  graph  is  geophysical  parameter  space  plot  from  the  principal  profiles 
(center  plots)  across  the  ordnance 
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y(m) 


105  mm  Artillery  Projectile 
Depth=0.49m,  !ndination=0,  Azimuth=90 


105  nm  Artillery  Projectile 
Depth=0.49m,  lndination=0,  Az]muth=0 


Figure  36.  Parameter  space  plots  for  TFM  and  IDEM  datasets  from  NRL  ordnance  signature  library. 

Top  two  sets  of  plots  are  for  105-mm  projectiles  with  identical  depths  and  inclinations  but 
different  azimuths  (0  and  90  deg).  Bottom  two  sets  of  plots  for  60-mm  (left)  and  81 -mm 
(right)  mortars  with  identical  depths,  azimuths,  and  inclinations 
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Figure  37.  NRL  MTADS  surveys  (TFM,  TFG,  TDEM)  over  cluster  of  eight  60-mm  mortars 


anomalies  caused  by  the  shallow  buried  60-mm  mortars  (Figure  38).  All  eight 
anomalies  caused  by  the  mortar  rounds  are  shown  dramatically  in  the  TDEM 
data,  while  two  of  Ae  eight  mortars  are  barely  visible  in  the  TFM  and  TFG  data 
with  the  selected  plotting  thresholds  and  contour  intervals.  The  TDEM  dataset 
is  used  to  estimate  location  and  depth  for  the  eight  anomalies.  Locations  are 
estimated  at  the  centers  of  the  TDEM  anomalies.  Depths  are  estimated  from 
anomaly  magnitudes  and  spatial  wavelengths,  knowing  that  all  targets  are 
60  mm  and  buried  at  shallow  depths.  Orientations  for  the  eight  targets  are 
estimated  from  characteristics  of  the  magnetic  anomalies,  and  for  illustrative 
purposes,  all  mortars  are  estimated  to  lie  at  an  inclination  of  0  deg  relative  to 
horizontal  and  to  have  azimuths  of  0,  45,  or  90  deg.  Using  the  target  estimates 
and  the  known  Earth’s  field,  magmod  is  used  to  compute  the  magnetic 
signatures  of  the  eight  mortars. 

To  compute  the  combined/composite  total  magnetic  field  of  the  eight  mor¬ 
tars,  a  MATLAB®  script  combine!  is  executed.  Each  case  is  computed  sepa¬ 
rately  using  magmod  and  each  saved  with  distinct  output  file  names.  The  eight 
output  files  from  magmod  are  input  to  combine!,  which  computes  the  composite 
field.  The  composite  field  output  from  combine!  is  subsequently  input  to  the 
magmod  to  display  the  output  composite  file.  After  minor  manual  iterations 
with  magmod  the  composite  predicted  TFM  field  is  shown  in  Figure  39.  The 
key  details,  of  the  TFM  map  of  the  60-mm  mortars  in  Figure  37,  are  replicated 
in  Figure  39. 
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Figure  39.  Result  of  manual  implementation  of  simplified  multisensor  integration  algorithm  (Figure  38)  to  replicate  total  magnetic  field  of  cluster  of  eight 
60-mm  mortars  measured  by  NRL  MTADS  system  at  Twentynine  Palms.  Ordnance  locations  and  depths  estimated  from  TDEM  images 
(Figure  37).  Initial  orientation  estimates  from  TFM  images,  then  manual  iteration  of  magmod  to  replicate  general  features  of  TFM  images 


7  Considerations  of 
Emerging  Technology 
for  UXO  Applications 


Introduction 


The  potential  importance  of  innovative,  emerging  geophysical  technologies 
for  advancing  UXO  discrimination  and  identification  or  classification  is  dis¬ 
cussed  in  Chapters  1  and  2  and  is  illustrated  in  this  chapter.  New  and  emerging 
technologies  that  may  allow  improved  discrimination  and  ultimately 
identification  of  UXO  are: 


a.  Multichannel,  multicomponent  TDEM  systems. 

b.  Multifrequency,  multicomponent  FDEM  systems. 

c.  Multicomponent  vector  magnetometers. 

d.  Tensor  magnetic  gradiometers. 

e.  Near  real-time,  high-resolution  microgravimetry. 

/.  Acoustic  wavefield  imaging. 

In  this  chapter,  the  potential  contribution  of  emerging  technologies  a  and  b 
above  is  assessed  by  examples  from  data  acquired  at  the  DARPA  sites.  The 
emerging  technology  TDEM  and  FDEM  systems  applied  at  the  DARPA  sites 
were  prototype  systems;  and,  as  a  result  of  the  lessons  learned  from  these  first 
data  collection  efforts,  needed  improvements  are  identified.  Some  data, 
particularly  for  the  prototype  TDEM,  were  noisy  and  do  not  represent  the  ulti¬ 
mate  potential  of  the  technology.  In  this  chapter,  selected  display  and  analyses 
approaches  are  explored  for  the  data  from  the  prototype  systems.  In  no  way  is 
the  total  potential  of  the  emerging  technologies  for  UXO  discrimination  and 
identification  assessed  or  achieved  in  this  report. 
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Multichannel,  Multicomponent 
IDEM  System 


Background 

A  prototype  multichannel,  multicomponent  TDEM  system,  designated  the 
EM61-3D‘,  was  used  for  data  acquisition  at  the  four  DARPA  sites  (Parsons 
Engineering  Science,  Inc.  1997).  The  EM61-3D  utilizes  the  standard  EM61  Tx 
coil  (1  m  X  1  m)  with  three  orthogonal,  0.5-m-diam  Rx  coils  within  the  Tx,  to 
measure  the  vertical  and  two  horizontal  components  of  the  secondary  induced 
transient  field  (i.e.,  a  multisensor  system).  The  system  uses  a  standard  TDEM 
receiver  electronics  package  which  samples  the  transient  with  20  channels 
(gates),  and  records  a  full  set  of  60  measurements  at  each  surface  location  (0.4-s 
intervals).  Figure  40  shows  the  EM61-D  in  operation  at  Fort  Carson, 

Colorado,  and  gives  the  times  of  the  20  gates  at  each  of  the  two  base  operating 
frequencies  (repetition  rate)  of  the  system.  Most  of  the  surveys  were  conducted 
at  the  7.5-Hz  repetition  rate  and  transients  (decay  curves)  from  three  Tx 
waveforms  were  stacked  to  produce  each  complete  record. 

Data  acquired  with  the  EM61-3D  over  one  of  the  1-hectare  DARPA  sites 
results  in  a  very  large  data  volume.  The  present  effort  concentrated  on  ways  of 
looking  at  small  segments  of  the  data  volume  to  assess  interpretive  procedures 
for  UXO  detection,  discrimination,  and  identification.  All  examples  presented 
below  are  from  surveys  at  the  Seabee  and  Turkey  Creek  sites  at  Fort  Carson, 
Colorado. 


Single-component  analyses  procedures 

Single-component  maps.  One  obvious  way  of  examining  the  data  is  to  pro¬ 
duce  site  maps  of  individual  components  for  selected  time  gates.  For  example, 
channels  (time-gates)  2  to  5  span  the  standard  EM6 1  window,  and  channel  3  is 
approximately  in  the  center  of  the  integrated  time-gate  recorded  by  the  standard 
EM6 1 .  Anomaly  features  are  observed  in  the  time-gate  maps,  including 
registration  targets  and  other  known  features.  There  is  obviously  background 
response  levels  due  to  the  geology  of  the  sites,  cultural  clutter,  and  a  significant 
level  of  system  noise;  this  suggests  applying  a  low-cut  amplitude  threshold  to 
enhance  anomalies  caused  by  buried  metallic  objects  (similar  to  the  CIS  layer 
query  in  Figure  17).  For  the  three  components  and  all  time-gates,  there  will  be 
60  site  maps  for  each  site.  With  this  rich  data  set,  there  are  considerable  possi¬ 
bilities  for  image  processing  type  analyses  for  background  subtraction  and 
localized  anomaly  enhancement,  change  detection  as  a  function  of  time,  and 
“false-color”  type  images  which  combine  the  three  components  at  selected 
times. 

Three-axis  representation  of  single-component  complete  decay  curves.  A 
second  procedure  is  to  examine  profiles  of  the  complete  time  decay  for  selected 
profile  lines.  The  three-axis  representation  in  Figure  41  is  for  the  vertical 
component  recorded  along  segment  of  a  north-south  line  at  the  Seabee  Site.  The 


^  Geonics,  Limited,  Mississauga,  Ontario,  Canada. 
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EM61-3D  Time  Gate  Locations  (in  ms) 
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Figure  40.  Photograph  of  the  EM61-3D  in  use  at  Fort  Carson,  Colorado  (top),  and  listing  of  the 

20  time  gates  for  each  of  two  operating  base  frequencies  (repetition  rates)  of  the  system 
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time  axis  “comes  out  of  the  plane  of  the  figure”  and  represents  approximately 
21.9  ms  of  the  decay  curve.  Noise  both  spatially  and  temporally  is  obvious  in 
the  representation,  however  anomalies  caused  by  known  buried  items  are 
apparent  and  two  of  the  knovm  (and  openly  releasable)  ordnance  items  are 
indicated.  Three  facts  are  evident  from  examination  of  numerous  raw  and 
spatially  and  temporally  smoothed  representations  like  Figure  41; 

a.  Buried  metallic  objects  have  higher  initial  (early-time)  TDEM  magnitudes 
than  background. 

b.  The  decay  persists  with  higher  magnitudes  to  later  times  than  background. 

c.  The  spatial  extent  of  TDEM  anomalies  along  track  are  greater  for 
ordnance  items  oriented  along  track  than  for  items  oriented  cross  track. 


Two-axis  representation  of  individual  component  complete  decay  curves. 
Another  technique  for  displaying  the  similar  information  to  that  in  Figure  41  is 
to  view  the  decay  curves  in  two  dimensions  as  illustrated  in  Figure  42.  In 
Figure  42  each  of  the  three  recorded  components  is  displayed  separately  in  a 
magnitude  versus  profile  distance  plot.  At  each  profile  position,  a  series  of 
vertical  dots  represents  the  20  measurements  of  magnitude  versus  time  of  the 
decay  curve.  The  background  clearly  forms  a  horizontal  “noise”  band;  the 
noise  band  is  considerably  smaller  for  the  y-component  (along-track  compo¬ 
nent).  Data  spikes  which  extend  above  or  below  the  noise  band  are  indicated  as 
targets  (T),  which  are  very  likely  buried  metallic  objects  and  possible  ordnance 
items.  Seven  targets  are  indicated  in  both  the  z-component  (vertical)  and 
y-component  data  at  the  same  locations,  while  three  of  the  same  targets  are  indi¬ 
cated  in  the  z-component  data.  Note  that  values  can  be  negative  for  the  x-  and 
y-components,  while  the  values  are  positive  for  the  z-component. 

Characteristics  of  complete  decay  curves  at  single  locations.  The  forms  of 
the  z-component  decay  curve  at  given  surface  locations  over  known  objects  are 
illustrated  in  Figure  31  (an  aluminum  plate)  and  in  Figure  43  (60-mm  and  81- 
mm  mortars,  105-mm  and  152-mm  projectiles).  As  suggested  by  McNeill  and 
Bosnar  (1996)  and  McNeill',  there  is  considerable  expectation  that  character¬ 
istics  of  the  complete  decay  curves  over  metallic  objects  will  provide  discrim¬ 
ination  capability  and  possible  identification  potential  for  UXO.  Several 
observations  supporting  the  preceding  proposition  are  apparent  in  Figures  31 
and  43  and  in  many  additional  decay  curves  examined  during  the  investigations 
(many  over  “bare”  ground  and  buried  non-metallic  objects): 


a.  All  of  the  metallic  objects  have  high  initial  (early)  time  magnitudes,  and 
after  an  initial  rapid  decay,  maintain  finite  magnitudes  for  late  times. 

b.  Nonmetallic  objects  have  low  to  intermediate  initial  magnitudes  and 
decay  to  near-zero  magnitudes  generally  in  less  than  5  ms. 

c.  The  decay  curve  for  the  aluminum  plate  decays  monotonically  to  very 
near-zero  magnitude  at  late  times. 


'  Personal  Communication,  July  1997,  J.D.  McNeill,  Geonics,  Lunited,  Mississauga,  Ontario, 
Canada. 
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Figure  41 .  Three-axis  portrayal  of  EM61-3D  vertical  component  raw  data  from  the  Seabee  Site,  Fort  Carson,  along  a  N-S  line.  The  designations  for  the 
orientations  of  the  two  identified  ordnance  items  are:  H— horizontal;  A-along-track;  and  C— cross-track 
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Figure  42.  Three-component  profile  representations  of  TDEM  decay  curves  at  Seabee  Site,  Fort  Carson;  20  decay  curve  measurements  plotted  at  each  profile 
position.  Targets  (T)  indicated  where  peaks  rise  above  the  “noise”  band  for  each  component 


Amplitude,  mv 


Figure  43.  IDEM  decay  curves  (vertical  component)  over  four  shallow  ordnance  items, 
Seabee  Site,  Fort  Carson 
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d.  The  decay  curves  for  the  ordnance  items  (ferrous)  are  quite  smooth  and 
indicate  nonmonotonic  decay  at  intermediate  times. 

e.  The  magnitude  at  channel  (time-gate)  20  increases  as  size  of  the  ordnance 
item  increases. 

The  above  observations  lend  considerable  support  to  the  potential  for  discrimi¬ 
nation  and  identification  of  UXO  from  observations  of  complete  decay  curves 
for  TDEM.  McNeill  and  Bosnar  (1996)  suggest  that  the  decay  curves  can  be 
fitted  with  various  forms  of  decaying  exponential  functions  with  two  to 
three  fitting  parameters,  which  may  be  diagnostic  for  ordnance  versus  non¬ 
ordnance  discrimination  and  of  ordnance  type. 

Single  parameter  discriminant.  Examination  of  Figures  31  and  43,  as  well 
as  many  other  decay  curves,  suggests  that  a  parameter  such  as  the  area  under  the 
curves  might  be  an  effective  discriminant.  The  area  is  related  to  the  “energy” 
induced  in  subsurface  volumes,  particularly  compact  metallic  objects.  Areas 
under  the  decay  curves  along  the  same  profile  line  shown  in  Figure  42  are 
displayed  in  profile  form  in  Figure  44.  The  same  seven  target  (T)  locations 
shown  in  Figure  42  are  indicated  in  Figure  44  and  correlate  with  dramatic  peaks 
in  decay  curve  areas,  with  the  background  areas  forming  a  very  narrow  “noise” 
band.  Although  not  examined  in  detail,  the  absolute  values  of  areas  of  x-  and 
y-component  decay  curves  above  buried  metallic  objects  apparently  will  also 
produce  maxima  above  a  background  noise  band. 


Multicomponent  parameter  space  analyses  techniques 

The  individual  three-component  displays  in  Figure  42  and  considerations  of 
the  decay  curve  areas  considerations  strongly  suggest  that  some  form  of 
multicomponent  (multisensor)  analysis  technique  will  provide  diagnostic  tools 
for  UXO  identification.  Two  types  of  multicomponent  parameter  space 
representations  are  presented  and  briefly  described  here:  (1)  spatial  profile 
figure  in  three-component  space  at  selected  times;  (2)  decay  time  figure  in 
three-component  space  at  selected  position.  The  first  type  multiparameter  space 
representation  is  a  spatial  profile  crossing  an  anomaly  feature  of  interest 
(identified  by  the  techniques  discussed  above),  where  each  “point”  on  the  figure 
represents  the  three-component  magnitudes  at  a  given  surface  location,  and  the 
entire  figure  is  for  a  selected  time  (i.e.,  there  could  be  20  such  figures).  The 
second  type  representation  is  for  a  selected  surface  location  and  the  figure  is 
formed  by  20  points,  with  each  point  representing  the  three-component 
magnitudes  at  a  specific  time.  Perhaps  the  greatest  difficulty  in  three- 
component  parameter  space  analyses  as  suggested  here  is  to  visualize, 
recognize,  and  characterize  the  geometry  of  the  figures  formed  by  the  two  types 
of  representations.  A  possible  approach  to  analyzing  the  geometry  is  to 
examine  projections  of  the  figures  on  the  three  orthogonal  component  space 
planes.  Another  procedure  for  geometric  and  time-evolution  analyses  requires 
reducing  the  space  to  2-D  by  combining  the  two  horizontal  components  into  a 
horizontal  magnitude,  as  was  done  for  Figtires  32  through  34. 

Spatial  profiles  in  three-component  space.  Examples  of  spatial  three- 
component  space  plots  are  shown  in  Figures  45  and  46.  Each  figure  shows 
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Figure  45.  Three-component  parameter  space  plot  for  a  background  location,  Seabee  Site,  Fort  Carson, 
for  EM61-3D  Channels  1, 10,  and  20 
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Figure  46.  Three-component  parameter  space  plot  for  an  aluminum  plate  (see  text  for  details),  Seabee  Site, 

Fort  Carson,  for  EM61-3D  Channels  1,10,  and  20 
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three  plots  for  Channels  1,  10,  and  20,  corresponding  to  early,  intermediate, 
and  late  times  (see  Figure  40  for  the  actual  times  of  the  channels).  Figure  45  is 
for  a  “background”  location,  i.e.,  a  location  with  no  known  buried  objects,  and 
hence  represents  the  response  of  4e  soil  at  the  site.  Tlie  axes  are  component 
magnimdes  in  mV,  with  the  z-component  the  vertical  axis.  The  15  points 
defining  the  geometric  figure  in  Figure  45  are  sequential  measurement  points 
along  a  N-S  profile.  Geometrically  the  figure  resembles  a  “scattergram”  that 
migrates  downward  and  to  the  left  as  a  function  of  time  in  the  plots. 

Ultimately,  at  very  late  time,  the  figure  should  shrink  to  a  point  at  (0,0,0);  the 
fact  that  it  doesn’t  appear  to  behave  in  this  manner  may  be  due  to  both  back¬ 
ground  (geologic)  variability  along  the  profile  line  and  system  noise. 
Examination  of  a  number  of  noisy,  fixed  location  decay  curves  (Figure  31  is  an 
example  of  a  nearly  noise-free  decay  curve,  wfrile  Figure  41  contains  some 
noisy  decay  curves) )  indicates  that  noise  in  the  vertical  component  may  be  as 
large  as  20  mV  for  some  channels  for  the  prototype  system.  Figure  45  indicates 
comparable  noise  levels  for  the  two  horizontal  components. 


Figure  46  shows  a  three-component,  profile  plot  crossing  a  20-  x  20-cm 
(8-  X  8-in.)  aluminum  plate  buried  at  a  depth  of  5  cm  to  the  top  of  the  plate. 
The  center  of  the  plate  is  crossed  near  profile  point  8.  Note  the  larger  axis 
intervals  in  Figure  46  than  in  Figure  45.  Geometrically,  the  plot  figure 
resembles  an  ellipse  with  a  cusp  at  the  origin,  similar  to  that  observed  for  some 
structures  in  gravity  gradient  space  plots  (Butler  1995).  The  plot  figure  shrinks 
in  size,  apparently  maintaining  the  cusp,  and  approaches  a  “point”  at  the  origin 
at  late  time  (Channel  20).  Actually,  considering  the  component  scale 
differences  between  Figures  45  and  46,  the  plot  figure  size  at  channel  20  is 
comparable  in  size  to  the  background  plot  at  Channel  20.  It  is  the  well-defined 
figure  shape  and  cusp  that  apparently  exists  for  the  aluminum  plate  that  holds 
promise  for  discrimination  and  identification  (Butler  1995). 


As  a  contrasting  case  to  the  conducting  flat  plate,  consider  an  iron  (ferrous) 
sphere,  i.e.,  magnetically  permeable  and  conducting  (Figure  47).  The  sphere  is 
12.4  cm  diameter,  buried  at  a  depth  of  10  cm  to  the  top,  and  lies  approximately 
beneath  profile  position  7  (indicated  by  an  asterisk).  An  important 
consideration  for  analyses  of  Figure  46  and  47  is  ttot  both  objects  are  buried  at 
depths  smaller  than  the  Tx  dimensions.  This  consideration  means  that,  except 
possibly  for  the  vertical  component,  symmetry  in  the  components  should  not  be 
expected,  as  it  would  for  a  deeply  buried  object,  vdiere  the  Tx  approximates  a 
dipole  source.  The  three-component  spatial  profile  plots  for  Channels  1,  10, 
and  20  for  this  case  (Figure  47)  indicate  a  much  different  geometrical  structure 
and  time  response  than  the  conducting  flat  plate  case.  From  the  perspective  view 
of  Figure  47  (identical  to  Figure  46),  the  plot  figure  appears  to  have  a  well- 
defined  geometry  at  early  time  but  Aen  collapses  apparently  to  a  “line”  at  late 
time.  The  point  directly  above  the  sphere  has  a  positive  cross-track  (x) 
component  at  early  time  that  becomes  negative  by  intermediate  times.  The 
cross-track  component  is  a  precise  indicator  of  the  location  of  the  object.  Noise 
is  clearly  contaminating  the  plot  figures  and  prevents  the  easy  recognition  of 
diagnostic  features,  particularly  from  the  one  perspective  view. 

Decay  times  in  three-component  space.  For  surface  positions  approx¬ 
imately  directly  over  selected  objects,  plate,  and  spheres,  time-decay  curves  in 
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three-component  space  are  plotted  in  Figure  48.  A  background  decay  curve  is 
also  shown  in  Figure  48  for  comparison.  The  aluminum  plate  and  iron  sphere 
are  the  same  targets  referenced  in  the  two  preceding  paragraphs.  The  aluminum 
sphere  (spherical  shell)  is  25.4  cm  in  diameter  and  is  bmied  60  cm  to  top. 

There  is  a  dramatic  difference  in  the  character  of  the  plate  and  sphere  responses. 
The  two  spheres  show  an  approximately  monotonic  decrease  (through 
intermediate  times),  in  the  vertical  (z)  and  cross-track  (x)  components,  and  then 
a  transition  back  to  increasing  values  of  the  cross-track  component.  Differences 
in  the  slopes  of  the  decay  curves  and  the  character  of  the  response  after  the 
transition  at  intermediate  times  are  caused  by  the  added  permeability  response 
for  the  iron  sphere.  Three-component  decay  plots  for  three  ordnance  items  (one 
mortar  and  two  artillery  projectiles)  and  one  metallic  AT  mine  are  shown  in 
Figure  49.  The  decay  responses  of  the  projectiles  are  totally  different  in 
character  from  the  decay  plots  in  Figure  48.  The  decay  response  of  the  metallic 
(steel)  AT  mine  is  totally  different  than  the  projectiles  responses  but  resembles 
the  aluminum  plate  response.  For  the  steel  mine,  there  is  also  a  transition  at 
intermediate  to  late  times  that  resembles  the  transition  of  the  iron  sphere. 


Multifrequency  FDEM  System 


The  second  emerging  technology  considered  as  part  of  this  effort  is  a  multi¬ 
frequency  FDEM  system  designated  the  GEM-3  (Geophex  1997).  The  system  is 
monostatic,  consisting  of  three  concentric  coils,  with  the  outer  coil  (nominally 
0.5-m  diam)  the  primary  Tx,  the  middle  coil  used  to  cancel  the  primary  field 
over  a  region  at  the  center,  creating  a  magnetic  cavity,  and  a  Rx  coil  within  the 
magnetic  cavity.  The  primary  field  is  generated  by  a  composite  waveform,  that 
can  generate  one  frequency  or  multiple  frequencies  simultaneously.  The  system 
is  broadband,  measuring  both  quadrature  and  in-phase  components  over  the 
frequency  range  90  Hz  to  24  kHz.  Data  acquisition  rates  up  to  10  Hz  are 
possible  with  multiple  frequencies,  and  the  system  can  also  operate  in  time 
domain. 


The  potential  of  the  multifrequency  FDEM  concept  is  illustrated  in  Fig¬ 
ures  50  and  51,  from  a  survey  of  the  Center  Square  area  of  the  Seabee  Site. 
The  ellipses  in  the  figures  surround  the  registration  targets  (an  aluminum  plate 
to  the  north  and  an  iron  sphere  to  the  south  at  each  location).  EMI  anomalies 
appear  in  the  plots,  for  4,050  Hz  and  12,270  Hz,  as  small,  tightly  grouped, 
closed  contours.  With  the  plotting  threshold  selected  for  the  anomaly  plots, 
most  background  noise  is  suppressed.  The  key  features  of  the  two  figures  are 
summarized: 


a.  For  the  anomaly  plot  of  the  quadrature  component  at  4,050  Hz  (Fig¬ 
ure  50),  only  the  southern  iron  sphere  targets  are  detected. 

b.  For  the  anomaly  plot  of  the  in-phase  component  at  12,270  Hz  (Fig¬ 
ure  51),  only  the  northern  aluminum  plates  are  detected. 
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Figure  48.  Three-component  time-decay  plots  for  background  (top  left),  aluminum  plate  (top  right),  iron  sphere  (bottom  right),  and  aluminxim  sphere 
(bottom  left).  Points  are  component  values  at  the  time  for  the  20  measurement  channels 
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Figure  49.  Three-component  time-decay  plots  for  8 1  -mm  mortar  (top  left;  inclination  90,  depth  to  top  20  cm),  1 05-mm  projectile  (top  right;  inclination  0, 

azimuth  90,  depth  to  top  55  cm),  152-mm  projectile  (bottom  left;  inclination  0,  azimuth  0,  depth  to  top  105  cm),  and  M15  metallic  (steel)  AT  mine 
(33.7  cm  diameter;  depth  to  top  15  cm).  jPoints  are  component  values  at  the  time  for  the  20  measurement  channels. 


,050  Hz,  Quadrature 
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Figure  50,  Multifrequency  FDEM  survey  results  at  Seabee  site.  Basic  site  layout  for  reference  (left)  and  an  EMI  response  map  (right)  at  4,050  Hz 
(quadrature  component)  for  a  portion  of  the  Center  Square  area  are  shown 


ase 
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Figure  5 1 .  Multi^frequency  FDEM  survey  results  at  Seabee  site.  EMI  response  maps  at  1 2,270  Hz,  in-phase  component  (left)  and  quadrature  component  (right). 
These  responses,  along  with  the  4,050ftz  response  in  Figure  50,  illustrate  material  type,  geometry  and  frequency  dependence 


c.  While  for  the  quadrature  component  at  12,270  Hz,  both  the  aluminum 
plate  and  the  iron  sphere  are  detected. 


These  features  indicate  that  the  EMI  response  is  both  frequency,  material  type, 
and  geometry  dependent. 

The  above  observations  suggest  the  possibility  of  electromagnetic  induction 
spectroscopy  (Private  communication,  I.  J.  Won,  Geophex,  Ltd.),*  where  the  in 
phase  and  quadrature  components  as  a  function  of  frequency  for  an  object  define 
a  EMI  spectra,  diagnostic  of  the  size,  geometry,  and  material  type  of  a  buried 
object.  Figure  52  illustrates  the  EMI  spectra  of  an  ordnance  item  (30  mm 
projectile)  and  various  metallic  objects  which  are  possible  clutter  sources  at 
sites.  Each  spectra  is  obtained  by  measuring  in  phase  and  quadrature 
components  at  a  single  location  above  the  objects  as  a  function  of  frequency, 
and  is  observed  to  be  distinctive  for  each  indicated  object.  Figure  53  carries  the 
spectra  observations  further,  by  indicating  spectra  dependence  on  orientation  for 
a  given  ordnance  item  (37  mm  projectile).  There  is  early  indication  that  the 
spectra  of  an  ordnance  item  with  inclination  I  ( 0  ^  I  ^  90  deg)  can  be  deduced 
from  in-plane  measurements  of  spectra  of  the  object  at  I  =  0  deg  and  I  = 

90  deg.  For  example,  the  spectra  for  the  case  I  =  45  deg  in  Figure  53  can  be 
replicated  by  simple  averaging  of  the  spectra  for  the  horizontal  and  vertical 
spectra  cases.  To  carry  this  concept  to  general  3-D  orientations  will  require 
three-component  measurements  and  development  of  data  transformations  from 
analytical  and/or  numerical  modeling  that  are  validated  by  measurements. 


Closure  for  Two 
Emerging  EMI  Technologies 


The  two  EMI  systems  briefly  examined  here  were  prototype  systems  at  the 
time  of  the  field  data  acquisition.  Both  systems  were  used  at  JPG  Phase  in 
UXO  Technology  Demonstrations  either  just  prior  to  or  just  after  their  use  at 
the  DARPA  sites,  and  both  systems  are  currently  in  varying  stages  of  advanced 
development  and  commercialization.  There  are  still  active  investigations  utiliz¬ 
ing  data  acquired  by  the  prototype  systems  at  the  DARPA  sites  and  at  JPG 
Phase  ni.  The  work  presented  here  and  also  the  work  of  other  investigators 
demonstrates  considerable  potential  for  UXO  discrimination  and  identification 
with  these  emerging  technology  systems.  The  full  potential  of  the  concepts  will 
likely  be  realized  only  if  the  systems  undergo  advanced  development  as 
multicomponent  and  multichannel  or  multifrequency  systems.  Also,  full 
phenomenological  modeling  capability  is  required  for  the  TDEM  and  FDEM 
cases.  In  spite  of  the  large  amounts  of  data  acquired  by  the  systems  compared 
to  standard  systems,  such  as  total  field  magnetometers,  the  standard  EM61,  and 
single-frequency  EMI  systems  and  metal  detectors,  both  systems  are  capable  of 
efficiently  surveying  large  areas. 


*  Geophex,  Ltd.,  Raleigh,  North  Carolina. 
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Figure  52,  Comparison  of  EMI  response  spectra  for  an  ordnance  item  (upper  left)  and  various  clutter  items; 
illustrates  discrimination  potential.  (Private  Communication,  1.  J.  Won,  Geophex,  Ltd.) 
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horizontal _  vertical,  nose  up 


Figure  53.  EMI  response  spectra  (multifrequency  FDEM)  for  a  37-mm  projectile  at  various  orientations 


8  Summary,  Conclusions, 
and  Recommendations 


Summary  and  Conclusions 


This  report  presents  a  broad  view  of  the  UXO  detection,  discrimination,  and 
identification  requirements  and  of  problems  associated  with  achieving  those 
requirements.  Concluding  that  UXO  detection,  except  for  problematic  cases,  is 
presently  possible  with  single-  or  multisensor  geophysical  datasets,  the  report 
concentrates  on  defining  the  current  status  of  detection  capability  and  then 
addresses  the  discrimination  and  identification  problems.  There  are  two 
problematic  cases  for  which  detection  capability  does  not  currently  exist: 

(1)  very  small,  widely  spaced,  single-ordnance  items,  such  as  20-mm 
projeailes;  (2)  ordnance  items  buried  too  deeply,  relative  to  their  size,  for 
detection  with  present  technology.  Very  small  ordnance  items  such  as  the 
frequently  present  20-mm  projectiles  are  generally  buried  at  very  shallow 
depths.  The  anomaly  magnitudes  and  spatial  wavelengths  of  such  small 
ordnance  items  prohibit  a  general  detection  capability  on  large,  live  sites  with 
current  technology,  although  closely  spaced  groups  or  clusters  of  small 
ordnance  items  may  be  detectable.  For  any  size  ordnance  and  site-specific 
conditions,  there  are  burial  depth  limits  below  which  the  object  cannot  be 
detected  for  each  current  technology  (e.g.,  TFM,  EMI,  GPR),  and  hence  there 
is  a  burial  depth  for  a  given  object  below  which  it  cannot  be  detected  by  any 
current  technology.  The  maximum  burial  depths  for  detection  by  current 
technology  are  fundamental  physical  limitations,  not  to  be  improved  or  solved 
by  improved  sensors  for  current  technology  (e.g.,  a  buried  object  at  a  “real 
world”  live  site  that  is  not  detected  with  a  TFM  sensor  with  0.01  nT  sensitivity 
will  not  be  detected  by  an  improved  TFM  sensor  with  0.001  nT  sensitivity). 

The  report  concludes  that  there  is  not  sufficient  information  content  in  cur¬ 
rently  deployed  single-sensor  survey  data  to  allow  UXO  discrimination,  except 
in  very  special  cases.  General  discrimination  capability  requires  multisensor 
datasets  and  integrated  or  joint  interpretation.  The  Naval  Research  Laboratory 
MTADS  system,  for  example,  can  acquire  high  quality  multisensor  (TFM, 
TFG,  and  TDEM)  datasets  at  high  data  densities  over  large  areas  at  reasonable 
production  rates  (current  production  rate  estimates  are  15  acres  (6  hectares  or 
0.06  km^)  per  day  for  TFM/TFG  and  7  to  8  acres  per  day  for  TDEM). 
Currently,  multisensor  integration  (“data  fusion”)  generally  occurs  at  the  data 
processing  and  interpretation  stage,  not  the  data  acquisition  stage.  Also, 
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multisensor  integration  by  most  practitioners  consists  of  side-by-side 
comparison  and  correlation  of  processed  data  maps  or  images.  Model-based 
interpretation  procedures  for  discrimination  are  rapidly  advancing,  but  are 
limited  to  analyses  of  each  type  of  sensor  data  separately,  followed  by  the  side- 
by-side  comparison  procedure  and  empirical  correlation  between  the  model 
interpretations.  There  are  also  efforts  to  achieve  multisensor  integration  using 
neural  network,  “fiizzy  logic,”  and  expert  system  type  approaches.  These  are 
generally  nonphysics-based  approaches,  and  they  rely  on  large  sensor  signature 
databases;  successful  application  of  these  approaches  for  UXO  discrimination 
are  reported.  This  report  assesses  procedures  for  model-based  interpretation 
and  presents  a  proposed  multisensor  data  integration  algorithm. 

The  report  also  concludes  that  UXO  identification  or  classification  is  not 
feasible  with  currently  deployed  technology.  Limited,  reported  successes  in 
UXO  identification  using  empirical  correlation  of  physical  properties  such  as 
interpreted  magnetic  moment  from  test  site  surveys  to  measured  magnetic 
moments  for  actual  ordnance  items  cannot  be  extended  to  general  surveys  of 
real  world  UXO  contaminated  sites.  A  key  component  of  the  ultimate  solution 
of  the  UXO  discrimination  and  identification  requirement  is  the  development  of 
physics-based  models  for  forward  and  inverse  modeling  of  geophysical  signa¬ 
tures.  A  distinct  possibility  exists  that  the  data  required  for  UXO  identification 
cannot  be  obtained  during  large  area  surveys  that  are  designed  to  detect  and 
map  UXO-like  targets.  The  logical  procedure,  following  detection,  is  to 
conduct  localized  area  surveys  over  the  identified  targets  (or  selected  targets) 
for  final  discrimination  and  identification. 

A  versatile  magnetic  modeling  program  for  UXO  is  documented  that  is  suc¬ 
cessfully  validated  by  comparison  to  measured  TFM  data.  The  program  models 
the  UXO  as  prolate  spheroids;  a  geometry  which  closely  approximates  the 
geometry  of  many  UXO  and  has  the  elongated  symmetry  axis  of  UXO.  The 
problem  of  TDEM  modeling  of  realistic  UXO  geometries  (cylinders  and  prolate 
spheroids)  is  considered  and  some  success  illustrated  for  the  full  time-domain 
solution.  However,  it  is  concluded  that  the  complete  time-domain  EMI  solution 
for  UXO  geometries  requires  an  extensive  and  dedicated  effort. 

Finally,  two  emerging  EMI  technologies,  a  multicomponent,  multichannel 
TDEM  system  and  a  multifrequency  FDEM  system,  are  evaluated.  Both  of 
these  systems  were  prototype  systems  at  the  time  of  acquisition  of  the  datasets 
used  for  this  study.  Various  analysis  techniques  for  the  data  from  these  proto¬ 
type  systems  is  considered.  The  potential  for  discrimination  and  identification 
with  these  systems  is  identified  and  is  considered  convincing  enough  to  recom¬ 
mend  continued  study  of  analysis  techniques  and  continued,  advanced  devel¬ 
opment  of  the  systems. 


Recommendations 


The  following  specific  recommendations  for  future  research,  development, 
and  technology  investments  are  the  mainly  the  result  of  the  work  documented  in 
this  report.  Some  of  the  recommendations,  however,  are  intuitive  and  already 
recognized  by  many  investigators. 
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Current  technology  multisensor  platforms.  For  cost-  and  time-efficient 
UXO  detection  surveys  of  large  areas,  a  self-propelled,  towed,  or  autonomous 
multisensor  platform  is  required.  The  platform  must  have  integrated  positioning 
capability.  The  platform  should  have  closely  spaced  arrays  of  sensors.  The 
recommended  sensors  for  the  platform  are  total  field  magnetometer  (TFM)  and 
time-domain  electromagnetic  induction  (TDEM)  systems.  The  TFM  arrays  can 
be  arranged  to  acquire  both  TFM  and  total  field  vertical  gradient  (TFG)  ^ta. 
Current  technology  may  require  separate  surveys  of  areas  with  the 
magnetometer  arrays  and  the  TDEM  arrays.  The  Navy’s  MTADS  system  is  an 
example  of  this  type  recommended  system  using  current  deployable  technology. 
Hand-carried  versions  of  the  current  technology  sensors  must  be  available  for 
use  in  densely  wooded  areas. 

Enhancements  for  current  technology  multisensor  platforms.  Readily 
identifiable  technology  enhancements  for  current  technology  multisensor  plat¬ 
forms  include;  (a)  greater  numbers  of  closer-spaced  sensors  in  the  arrays, 
particularly  to  address  the  very  small-size  ordnance  problem;  (b)  sensor  systems 
with  greater  sampling  rates  along  track;  (c)  capability  to  acquire  TFM  and 
TDEM  data  simultaneously;  (d)  “real-time”  UXO  detection  decision  aids.  The 
development  and  utilization  of  UXO  test  sites  is  a  necessary  adjunct  to  the 
recommended  enhancements. 

Multisensor  data  integration  for  UXO  discrimination.  Except  for  posi¬ 
tioning  data,  all  integration  should  be  accomplished  at  the  data  processing  and 
interpretation  stage  and  not  at  the  data  acquisition  stage.  While  there  may  be 
some  potential  in  the  approach  of  using  neural  networks  and  other  nonphysics- 
based  approaches  for  multisensor  integration  for  discrimination,  the  recommen¬ 
dation  from  the  present  work  is  that  data  integration  and  interpretation  be 
achieved  by  fundamental  physics-based  modeling  approaches.  The  recommen¬ 
dation  does  not  preclude  the  utilization  or  fundamental  value  of  large,  measured 
UXO  signature  libraries  for  validation  of  modeling  capability,  or  the  use  of 
large  libraries  of  model  computed  UXO  signatures  for  signature  matching 
approaches  to  discrimination  and  identification. 

Emerging  technology  for  UXO  discrimination  and  identification.  This 
report  identifies  several  emerging  technologies  in  Chapter  7  and  discusses  two 
of  the  technologies  in  detail.  Emerging  technologies  offer  the  greatest  hope  and 
potential  for  capability  for  remote  discrimination  of  UXO  anomalies  from  false 
alarm  anomalies  and  for  identification  or  classification  of  buried  UXO. 

Research  and  development  of  these  emerging  technologies  must  be  aggressive, 
with  the  goal  of  integrating  the  new  capabilities  into  deployable  systems  and 
platforms. 

Integration  of  new  technology  into  deployable  systems  and  multisensor 
platforms.  Integration  of  new  technological  capability  into  deployable  hand¬ 
held  systems  and  multisensor  platforms  must  be  encouraged  and  supported,  but 
only  after  needed  research  and  development  and  field  validation  at  well- 
controlled  UXO  test  sites.  Accompanying  the  emerging  technology  sensor 
development  must  be  data  processing  and  interpretation  procedures  to  complete 
the  technology.  Some  of  the  emerging  technologies  may  prove  to  have  general 
applicability,  while  some  will  be  limited  to  localized  investigation  for  UXO 
identification.  Also,  some  may  prove  to  have  very  limited  applicability  and  not 
justify  any  type  advanced  development.  Some  of  the  identified  emerging 
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technology  may  be  available  for  integration  into  deployable  systems  in  the  near- 
term,  say  2  to  3  years:  e.g.,  (a)  three-axis  vector  magnetometers,  (b)  multiple 
axis  magnetic  gradiometers;  (c)  single-component  (vertical)  multichannel 
TDEM  systems;  (d)  single-component  multiffequency  FDEM  systems;  (e)  near 
real-time,  high-resolution  microgravimetry.  Other  emerging  technologies  may 
require  a  longer  term  investment,  say  3  to  5  years,  such  as  multicomponent 
TDEM  and  FDEM  systems  and  &11  tensor  magnetic  gradiometers.  The  final 
identified  emerging  technology,  acoustic  wavefield  imaging,  may  not  prove 
feasible  or  applicable  to  real  world  UXO  problems  and  certainly  represents  a 
longer  term  investment. 
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Data  cor  to  1  B.S.:60mm,  60  cm  offset,  all  orientations 
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pb4=uicontrol(main,'Style','push',. . . 

'Position', [40,260, 245,40]*sss,...  'load  lastp.mat;namep=namef;',... 

'String','PLOT  RESULTS',...  'save  lastp.mat  dipole  octopole  dattyp  plttyp  namep  oval;',. 

'CallBack',['load  lastp.mat;figure(h);',. . .  'set(edxx,"string",namep);figure(main);']); 
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'magfield(namef,Xmax(  1  ),Xmin(  1 ),  Ymax(l ),  Yniin(  1  ),Z(  1  ),B0mag(  1  ),ldeg(  'BackgroundColor','b'); 

l),Ddeg(l),',...  ; 

pb_close=uicontrol(infig,'Style','push',. . 

'L(l),D(l),div,perm(l),Thetadeg(l),Phideg(l),XO(l),YO(l),ZO(l),MorF);',..  'Position',[20,540,100,30]’»sss,... 
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'Representation  of  Spheroidal  Body|',. . .  yy=yy-80;dostuff; 

'Surface  of  Spheroidal  Bodyj',...  ndogs=4; 

'Listing  of  Major  Parameters'];  strl  a=num2str(cval(l  ,ndogs));str2a=num2str(cval(2,ndogs)); 

str3a=num2str(cval(3,ndogs));str4a=nuin2str(cval(4,ndogs)); 
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'],...  'UserData'.ncat,... 

'CallBack',['mode=get(typmen,"Value");']);  'BackgroundColor'.'w',... 

;  'CallBack',['forkcat=l:16,'.... 

txq=uicontrol(outfig,'Style','textVUnitsVpixels',...  'vvv=get(ed(kcat), "String");',. 

'Position'.[460,540,130,20]*sss,...  'vvvv=([ "["  ,vw, "]" 


'ifkcat=l,vl=str2num(vvw);end;',...  measfile=compfile; 

'if  kcat“2,v2=str2num(vvw);end;',. . .  readmeas=0; 

'if  kcat— 3,v3=str2num(ww);end;',...  if  exist(measfile)~2,readmeas=l  ;end; 

'if  kcat=4,v4=str2num(vvw);end;',. . .  ; 

'ifkcat=5.v5=str2num(vvwy,end;’,...  OUTX=zeros(nx,ny);OUTY=zeros(nx,ny);OUTZ=zeros(nx,ny); 


i-T  5^  .  ^  ^  ^  m  rf  m  NO  00  ^ 

E  ^  >  CO  .>  o  O  r-  <s  -jj  -j|  -||  11  ll  1l  ll  II  II  II  ll 


'O'  vO  w'  'w' 


•  -  •  •  TS  T3  'TI3  *T3 

d  C5  cS  C  c 
•o  *0  -b  -b  «  <J>  H  «  «  « 

C  C  C  C  /i-N 

<0^0  ^^  >>>>>> 


C  G  z:  G  G  G  q 

^  ^  ^  y  <2 

b  ^  b  ^  ^  <b  ^ 

Vi  tn  tn  In  V3  ui  tn 

II  II  II  II  II  11  II 

o  ^  «S  in  VO 

>. 

,  o  ^  <s  fO  m  VO 


03  03  A  OS  a  eo  a 

o  o  o  o  o  o  o 

5tab3?tab:b35ta5ta-b 

. § 


IS:'*-- 

S  o  o  '5' 

i§,|  I 

is  i-s 

Q  55  g  -q 

g  O  ® 

OQ  00 


It 

lib- 

Cu  S 
00  P 

41 

ft"  a 


Appendix  B  MAGMOD  Program  Listing 


B7 


o 

H 

O 

Cu, 


I 

i! 

§)  § 

Cb  S 

I S 
s  §■ 

•3  I 


W 

<u 

1 

o. 

g 


o 

g 


8 
S' 

Cu 

U  W3  ^ 
N©  nP  ^ 

^  ^  t 

x> 
'B 

(tT 

CO 

P 


J  a 

II  I 


•r-Sa  g 

bO  o  o 

yO  V©  <vO  %©  <<xP  >,^  s,©  >yO  N,©  «sP 

©^  O*^  ©^  ©^  ©^  ©^  ©^  ©^  ©^  ©''  ©^  O  0'‘ 


-73 

§ 

f 

Vi 

B 

*0, 

c 

\o 

4  - 

^  !l 
g- 1 

V  ?t3  &■ 


O 

IM  4> 

<2 

If 
i  ^ 

4> 

c 

c: 

4> 

Q 


.  ^  -O  T3 

^2  cj  c 

c  flj  a> 

4)  i*V  ^ 
•  -  CS 

2  %  ^ 
8  I  I 

CT  O^  CT 

^  cs 

AAA 

.  Vi  tn  Vi 

V  S  S  S 
'o-’cL’H. 
S  c  c  c 
a^5ta  5t2  hd 


«Wh 

sp 

-o" 


^  1 


lEi: 

Is 

B  ^ 

b  II 


i,^ 

•'  IT' 

il 

■c  s 

.1^ 


■s 

a. 

Q 

H  Vi 

i_j  ^ 

il  O 


Qu  4J 

<U  U4 

Q 


|... 
ifil  •. 

lill 
1 1  i  i  . 
gllff 
«g  8  I  S 

jr  V  g  o  o 

S  £  g  c  ^  -g 

<l> 


•o 

g 


•i  ^ 
s  s 

TO  03 


1l  1l 


ffi 


^  CS  CO 
CS  CM  CS 

'w'  %_✓ 

11^ 

*0  "Td  "TD 


II  II  >«,,' 

/— s  O 

oo  o\  ^ 


"  “><  isj  G  ia  > 

1l  1)  1l  1l  1l  1l  11 

y-N  y-v  y— s  y— > 

^  cN  CO  \o 

’■^-y  v-y  v-y  s— 'w-y  's-y  v»y  -w 
CM<NCN<NCNCM(N<N<N<N 

1111111111 

.  «^  ^O  *0  *t3  *tD  T3  *13  *T3  "O  "O  *0 


!!, 


_o 


II 

o 


5 

« 

•o 

I 

o 

Bu 
.  c 


«3 

CO 

s 

2. 

*5 

> 

Vi 

Vi 

2 

U 

3 

jh 

"o 


JL’g'S 

-y  p 

II J  3 

^  ^  ^  -O  i 

?llt|| 

i  o  'tj 


& 


.g 


B8 


Appendix  B  MAGMOD  Program  Listing 


ip=l;  mz:3:nx,mniz:3:ny)), 

bodyproj;  %index  for  profile  plot  matrix.  hold  on;contoiir(X,Y,OUTanom);for  bbbl=l  :lxO, 

for  mp=l  :nplots,  %Step  through  the  subwindows.  plot(Xrot(bbbl,:),Yrot(bbbl,:),'W:');end;holdofif; 

if  datatype(mp)=-b'  |  datatype(mp)— 'B',  title('Horizontal  Vectors');axis('square');axis('equar); 

colormap('default');  xlabel(Northlb);ylabel(Westlb); 
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B9 


elseifdatatype(mp)=-v';  xlabel(Northlb);ylabel(Westlb);zlabel('nanoteslas'); 

subplotx(mode,qrow,qcol,mp),bodyproj ;  elseif  plottype(mp)=='M' 

minz=rem(nfi-ame,3)+l ;  subplotx(mode,qrow,qcol,mp),meshc(X,Y,dfile), 

quiver(X(mmz:3  :nx),  Y(rmnz:3  :ny),0UTX(mmz:3  :nx,mmz:3  :ny),OUT  Y  (m  ; 


axax=axis;  ifnframe===K),nix=TOax(max(dfile));nm=min(inin(dfile)); 

maxmax==max(max(dfile));minmin=min(min(dfile));  vcx)ntour=[nin:(mx-mn)/20:nix]; 

text(Xniax,Yinin,axax(6),['  Max  ',num2str(maxmax),'  Min 

’,num2str(minmin)]);  conlabl=contour(X,Y,dfile,vcontour);axis(’square');axis('equar);%clabel(con 
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title^label);xlabel(Westlb);ylabel('nanoteslas');  ax=axis;range=ax(4)-ax(3);axis([ax(l),ax(2),ax(4)-2*range,ax(4)]); 

ax=axis;range=ax(4)-ax(3);axis([ax(  1  ),ax(2),ax(4)-2*range,ax(4)]);  hold  on; 

hold  on;  if  readmeas==  1 

if  readmeas=l  plot(measured(: ,  1  ),measured(:,2),'w+'); 
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holdofif;  title(plabel);xlabel(Westlb);ylabel(’nanoteslas'); 

ip=ip+l  ;  hold  on; 

elseif  plottype(mp)='$';  if  readmeas==l 

subplotx(mode,qrow,qcol,mp);  plot(measured(:,l),measured(:,2),'w+’); 
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%  center  location.  %  This  MATLAB  script  calculates  the  dipole  and  octopole 

[xsph,ysph,zsph]=sphere;  %Create  a  default  sphere.  %  contributions  to  the  magnetic  field  at  each  grid 

xsph=D(bbb  1  )*xsph/2;  ®>^Scale  the  x  axis  of  spheroid.  %  point  for  the  body  and  field  data  specified  in  SPHEROID .M 

ysph=D^bb  1  )*ysph/2 ;  %Scale  the  y  axis  of  spheroid.  % 
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DIPy=zeros(nx,ny);  %  back  to  surface  coordinates  as  Bd. 

DIPz=zeros(nx,ny);  Bd=Rot'*bd; 

DIPmag=zeros(nx,ny);  DIPx(n,m)=Bd(l);  %fill  in  X  component  of  dipole  field. 

DIPanom=zeros(nx,ny);  DIPy(n,m)=Bd(2);  %fill  in  Y  component  of  dipole  field. 

OCTx=zeros(nx,ny);  DIPz(n,m)=Bd(3);  %fill  in  Z  component  of  dipole  field. 
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)*termB/rmagM)/(4*pi*rmag''5);  OCTz=OCTz';  %  orientation  for  plots. 

%end  octopole  calculation  OCTmag=OCTmag'; 

%  OCTanom=OCTanom'; 

%  Transform  magnetic  dipole  field  bd  from  body  coordinates  eval(['save  '.name,'  MorF  BOmag  Ideg  Theta  Thetadeg  Phi  Phideg  L  D  perm 
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Z1=Z;  Y0s=[Y0s  YO]; 

Z0s=[Z0sZ0]; 

%  Zero  the  values  to  be  summed  Ds=|T)s  D]; 

;  Ls=[Ls  L]; 


Phidegs=[Phidegs  Phideg];  DIPx=DIPxs; 

Thetadegs=[Thetadegs  Thetadeg];  DIP)F'DIPys; 

perms=[penns  perm];  DIPz=DIPzs; 

OCTanom=OCTanoms; 

DIPanoms=DIPanoms+DIPanom;  OCTmag=OCTmags; 
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fprintf(ihand,'\nV  INPUT  PARAMETERSW);  distance82=sqrt((X(lmxx82)-X(linnx82)r2-KY(linxx82)-Y(lmny82))''2); 

fprintf(ihand,['  File  ',namein,'\n\r  Date  ^day.'NnVr  Time'.t,' \n\r  \n\r']);  ; 

Q)rintf(ihand,['  UXO  Location(x,y,z)  %10.4f%10.4f  lprmtf(ihand,' RESULTS  Total  Field  AnomalyVnV); 

%10.4f\n\r'],X0,Y0^);  $rintf(ihand,'  Max  X  Y  Min  X  Y  DistW); 

fprintf(ihand.['  UXO  Size(length,diameter)  %10.4f%10.41\nV],L,D);  ^rintf(ihand,['  OCT%10.2f%6.2f%6.2f%10.2f%6.2f%6.2f 
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[maxoctdip,lnixy82]=max(max(OCTDPanom'));  ^rintf(ihand,'  RESULTS  Magnitude\n\r'); 

[minoctdip,lnmy82]=min(min(OCTDIPanom'));  Q)rintf(ihand,'  Max  X  Y  Min  X  Y  Dist  \n\r'); 

distance8=sqrt((X(lmxx8)-X(lninx8)r2+(Y(lmxx8)-Y(linny8))''2);  4)rintf(ihand,['  OCT%10.2f%6.2f%6.2f%10.2f%6.2f%6.2f 

distance2=sqrt((X(lmxx2)-X(linnx2))''2+(Y(ltnxx2)-Y(linny2)r2);  %9.21\nV],.. 
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[maxdip,lmxx2]==max(inax(DIPmag));  [maxcx:tdip,lmxx82]=niax(max(OCTDIPx)); 

[mindip,lmnx2]==miin(niin(DIPmag));  [minoctdip,lrnnx82]==niin(min(OCTDIPx)); 

[maxdip,Inixy2]=max(max(DIPmag’));  [maxoctdip,lmxy82]=max(max(OCTDIPx')); 

[mindip,lmny2]=min(minG!)IPmag’));  [mmoctdip,lniny82]=min(min(OCTDIPx')); 


distance8=sqrt((X(ltnxx8)-X(linnx8))''2+(Y(lmxx8)-Y(lnmy8))''2);  fprintf(ihand,'  Y  Component  Max  Min  Xloc  Yloc  Dist 

distance2=sqrt((X(lnixx2)-X(lmnx2))''2+(Y(lnixx2)-Y(lnmy2))''2);  \nV’); 

distance82=sqrt((X(lnixx82)-X(lmnx82))''2+(Y(lmxx82)-Y(lmny82))^2);  fprintf(ihand,['  OCT  %10.4f  %10.4f%10.4f\nV],... 

;  maxoct,X(lnixx8),Y(lnixy8)); 

Q)rintf(ihand,' X  Component  Max  Min  Xloc  Yloc  Dist  Q)rintf(ihand,['  OCT  %10.4f%10.4f%10.4f%10.4f\n\r'],... 
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distance8=sqrt((X(lmxx8)-X(lmnx8))''2+(Y(lmxx8)-Y(lmny8))^2);  fprintf(ihand,' Z  Component  Max  Min  Xloc  Yloc  Dist 

distance2=sqrt((X(lmxx2)-X(lmnx2))''2+(Y(lmxx2)-Y(lmny2))''2);  \n\r'); 

distance82=sqrt((X(lmxx82)-X(lmnx82))''2+(Y(lmxx82)-Y(lmny82))''2);  Q)rintf(ihand,['  OCT  %10.4f  %10.4f%10.4f\nV],... 

;  maxoct,X(lnixx8),Y(lnBgf8)); 
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iil5=length(Z0); 


[vm,lm]==max(dfile(cv(ip,indx),:)); 
[vn,ln]=min(dfile(cv(ip,indx),:)); 
difl(indx)=abs(  Y(lin)- Y  (In)); 
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[vm,lm]=max(dfile(:  ,cv(ip,indx)));  'String'.thedata,' Value', dattyp(ndogs),. . . 

[vn,ln]=min(dfile(:,cv(ip,indx)));  'CallBack',['for  kdogs=l  :6,dattyp(kdogs)=get(cb(kdogs),"Value");',. 

difl(indx)=abs(Y(lm)-Y(ln));  'if  dattyp(kdogs)=l  |  dattyp(kdogs)  >=  7,',... 

else  'set(pm(kdogs),"Visible","off ');',. . . 
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evl(ndogs)=uicontrol(outfig, 'Style', ’edit',...  else,set(ev2(ndogs),'Visible','off);end; 

'Position', [xxxx,yy,l  10,15]*sss,...  if  tval  >=  2,  set(ev3(ndogs),'Visible',’or 

'String',strla,'BackgroundColor','c',...  else,set(ev3(ndogs),'Visible','oflP);end; 

'Visible','off,...  if  tval  >=  3,  set(ev4(ndogs),' Visible', 'or 

'CallBack',['forkdogs=l:4,strla=get(evl(kdogs), "String");',...  else,set(ev4(ndogs),'Visible','ofir);end; 


if  mode=2 
xlenl=70; 
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yst2=[.2..60];  ystl=25; 

xst3=[.25,.6..25,.6];  xst2=[.  1 5..50]; 

yst3=[.2,.2,.60,60];  yst2=[.25,.25]; 

xst5=[.25,.6,.26,.6..26..6];  xst3=[.  1 5,.6,.  1 5,.6]; 


yst3=[.25,.25,.60,.60];  iil2=length(Phideg); 

xst5=[.  1 5,.  1 5..40,.40,.65,.65];  ii  1 3=length(X0); 

yst5=[.25,.6,.25..6,.25,.6];  iil4=length(Y0); 

end;  iil5=length(Z0); 
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ii7=length(Ideg);  ); 

ii8=length(L);  figure(h); 

ii9=length(D);  eval(['print  -dgifS  F',mim2str(iframe),'.GIF']); 

ii  1 0=length(perm);  rectzzz=get(gcf, 'Position'); 

iil  l=length(Thetadeg);  rectzzz(l  :2)=[0  0]; 


ififrame==0,  kk=2; 

clear  M;  ; 

M=zeros(size(getframe),numf);  while  kk==2, 

end; 

iframe=iframe+l ;  name=['m’,int2str(jj),’.inat'] ; 
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figure(main);  % 

q=7./6.; 

-»  yup=linspace(0, 1,10); 

xover=linspace(0,l,10); 

;  xorig=0*xover; 

jj=0;  yorig=0*yup; 
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if  octopole==l 

text(0.4,l  .8*q/This  is  only  the  octopole  contribution'); 
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